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1. Zur Quantentheorie der Molekeln; 
von M. Born und R. Oppenheimer 


Es wird gezeigt, daß die bekannten Anteile der Terme einer Molekel, 
die der Energie der Elektronenbewegung, der Kernschwingungen und 
der Rotationen entsprechen, systematisch als die Glieder einer Potenz- 
entwicklung nach der vierten Wurzel des Verhältnisses Elektronenmasse 
zu (mittlerer) Kernmasse gewonnen werden können. Das Verfahren 
liefert u. a. eine Gleichung für die Rotationen, die eine Verallgemeine- 
rung des Ansatzes von Kramers und Pauli (Kreisel mit eingebautem 
Schwungrad) darstellt. Ferner ergibt sich eine Rechtfertigung der von 
Franck und Condon angestellten Betrachtungen über die Intensität 
von Bandenlinien. Die Verhältnisse werden am Beispiel der zwei- 
atomigen Molekeln erläutert. 


Einleitung 

Die Terme der Molekelspektren setzen sich bekanntlich 
aus Anteilen verschiedener Größenordnung zusammen; der 
größte Beitrag rührt von der Elektronenbewegung um die 
Kerne her, dann folgt ein Beitrag der Kernschwingungen, 
endlich die von den Kernrotationen erzeugten Anteile. Der 
Grund für die Möglichkeit einer solchen Ordnung liegt offen- 
sichtlich in der Größe der Masse der Kerne, verglichen mit 
der der Elektronen. Vom Standpunkte der älteren Quanten- 
theorie, die die stationären Zustände mit Hilfe der klassischen 
Mechanik berechnet, ist dieser Gedanke von Born und 
Heisenberg’) durchgeführt worden; es wurde gezeigt, daß 
die aufgezählten Energieanteile als die Glieder wachsender 


Ordnung hinsichtlich des Verhältnisses ]/ “7 erscheinen, wo m 


die Elektronenmasse, M eine mittlere Kernmasse ist, Dabei 
traten aber Kernschwingungen und Rotationen in der gleichen 
(zweiten) Ordnung auf, was dem empirischen Befund (bei kleinen 
Rotationsquantenzahlen) widerspricht. 


1) M. Born u. W. Heisenberg, Ann. d. Phys. 74. S. 1. 1924. 
Annalen der Physik. IV. Folge. 84. 30 
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Hier soll die Frage vom Standpunkte der Quantenmechanik 
neu aufgenommen werden.') Dabei stellt sich heraus, daß 


man nicht nach V I sondern nach % a entwickeln muß, 
um die natürliche Reihenfolge der Energieanteile zu bekommen, 
Die Überlegungen werden überdies viel einfacher und durch- 
sichtiger als die der älteren Theorie. Die Kernschwingungen 
entsprechen den Gliedern 2. Ordnung, die Rotationen denen 
4. Ordnung in der Entwicklung der Energie, während die 
Glieder 1. und 3. Ordnung verschwinden. Das Ausfallen der 
Glieder 1. Ordnung hängt mit der Existenz einer „Gleich- 
gewichtslage“ der Kerne zusammen, in der die Elektronen- 
energie bei ruhenden Kernen ein Minimum ist. Die Terme 
4. Ordnung für die Rotationsbewegungen stellen eine Verall- 
gemeinerung des bekannten Ansatzes von Kramers und Pauli?) 
dar, die das Verhalten einer Molekel durch das eines Kreisels 
mit eingebautem Schwungrad wiedergegeben haben. Um die 
Eigenfunktionen und damit die Übergangswahrscheinlichkeiten 
auch nur in nullter Näherung zu bestimmen, muß die Energie- 
berechnung gerade bis zu den Gliedern 4. Ordnung (Rotationen) 
getrieben werden. Man erhält für die Wahrscheinlichkeiten 
gleichzeitiger Sprünge von Elektronen-, Schwingungs- und 
Rotationsquantenzahlen Ausdrücke, durch welche die von 
Franck’) entwickelten, von Condon‘) ausgestalteten Vor- 
stellungen eine präzise Fassung erhalten. 

Die Näherungen von höherer als 4. Ordnung werden in 
der vorliegenden Arbeit nicht behandelt; sie würden den Koppe- 
lungen zwischen den drei Grundtypen der Bewegung ent- 
sprechen. Eine Durchrechnung dieser Effekte hätte nur Sinn 
bei gleichzeitiger Berücksichtigung aller Entartungen der 
Elektronenbewegung bei ruhenden Kernen, vor allem der 
Heisenbergschen Resonanzentartung infolge der Gleichheit 


1) Bei der Diskussion der Grundlagen dieser Arbeit hat uns Hr. 
Dr. P. Jordan durch wertvolle Bemerkungen unterstützt, wofür wir 
ihm unsern Dank zum Ausdruck bringen möchten. 

2) H. A. Kramers, Zeitschr. f. Phys. 13. S. 343. 1923; H. A. 
Kramers u. W. Pauli jr. 13. S. 351. 1928. 

8) J. Franck, Trans. Faraday Soc. 1925. 

4) E. Condon, Phys. Rev. 28. 8. 1182. 1926; Proc. Nat. Acad. 13. 
8. 462. 1927. 
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der Elektronen (gegebenenfalls auch einiger Kerne) und bei 
zweiatomigen Molekeln der Entartung der Eigenrotation um 
die Kernverbindungslinie; auf diese verwickelten Betrachtungen 
soll hier verzichtet werden. 

Als Beispiel werden die zweiatomigen Molekeln näher 
betrachtet, und zwar nicht nur nach der allgemeinen Methode, 
sondern noch nach einer andern, bei der die Rotation durch 
Separation der Variabeln bereits in nullter Näherung berück- 
sichtigt wird, ähnlich, wie es Born und Hückel!) in der 
älteren Quantentheorie gemacht haben. 


§ 1. Bezeichnungen und Definitionen 


Die Masse und die rechtwinkligen Koordinaten der Elek- 
tronen bezeichnen wir mit 


Tey Yar Ze 
die der Kerne mit 
My X,, 4,. 
Ist M irgendein Mittelwert der M,, so setzen wir 


dann sind die u, reine Zahlen der Größenordnung 1. Die 
potentielle Energie des Systems sei 


= U(z, 2; 

hier, wie überall im folgenden, werden wir durch z die Ge- 
samtheit der Elektronenkoordinaten, durch X die der Kern- 
koordinaten darstellen. Die Funktion U hängt nur von der 
relativen Lage der Partikeln ab; doch machen wir von ihrer 
speziellen Formen (Coulombsches Gesetz) keinen Gebrauch. 
Der kinetischen Energie der Elektronen N der Operator 


1) M. Born u. E. Hückel, Phys. Ztschr. 24. S, 1. 1923. 
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wobei das Zeichen >; die Summe über Glieder bedeutet, die 


aus dem angeschriebenen durch zyklische Vertauschung vo 
x, y, z hervorgehen. 
Die kinetische Energie der Kerne ist 


h? 0° 
6) 
Die Gesamtenergie entspricht dem Operator 


(6) H=H,+»*H,, 
wo 


ae U= H,(2, 


U 
TE = x4 H, (sx) 
gesetzt ist. 
Wir führen nun statt der rechtwinkligen Kernkoordinaten 


neue ein, nämlich 3N—6 Funktionen 
(8) &= §,(X), 


die die relative Lage der Kerne zueinander festlegen, und 
6 Funktionen 


(9) = F(X), 

die die Lage der Kernkonfiguration im Raume bestimmen. 
Man kann das in symmetrischer Weise etwa so machen, dab 
man die rechtwinkligen Koordinaten r,, y,, 3, der Kerne be 
züglich der instantanen Hauptträgheitsachsen einführt; zwischen 
diesen bestehen dann die 6 Gleichungen 


=9,. 
Man kann die r,, ... also durch 3N—6 unabhängige Parameter 
&,, &,--. ausdrücken: 


5“ PAIR 
Sodann bestehen zwischen den ursprünglichen und den neuen 
Koordinaten Transformationsgleichungen der Form 


(10) X,= X, +34 ay", Ps W)Y, (§); 


X; 4 sind die und die «,, die 
Koeffizienten der orthogonalen Drehungsmatrix, also bekannte 
Funktionen der Eulerschen Winkel #, y, w. Die Größen 
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X: Yo, 2, +, 9, w sind die in (9) mit #, bezeichneten 
Funktionen. Durch (10) sind die X, als Funktionen der #,, &, 
bestimmt; durch Auflösen bekommt man die Ausdrücke (8) 
und (9). 

Durch diese Transformation wird die Energie 7 natürlich 
nicht in Anteile der Translation, Rotation und Relativbewegung 
der Kerne separiert. Wohl aber kann man H, in 3 Teile 
spalten, die sich im folgenden charakteristisch verschieden 
verhalten werden: 


(11) H, = H, + 4, 


H,, ist homogen linear in den 45 ; H,, enthält die 


i, ist von allen Ableitungen nach den &, frei. Man kann 
über diese Operatoren noch einige allgemeine Aussagen machen. 
Wendet man den ganzen Operator H, auf irgendeine Funktion 
f(§ der relativen Kernkoordinaten §, an, so muß die ent- 
stehende Größe H, f(£) von der Lage im Raume, also von den 
d,, unabhängig sein. Insbesondere können in H,. die Koeffi- 


zienten der nicht von den #, abhängen. Dagegen 
kommen in Are neben die &,, 9, und , 
neben den 38,55; die => > vor. 


Für sweiatomige Molekeln werden wir diese Operator- 
fanktionen explizite angeben. 
Das zu lésende mechanische Problem lautet 


(12) (H, + «4H, — W)y =0. 


Wir werden zeigen, daß diejenige Lösung, welche einer Ver- 
einigung der Kerne und Elektronen zu einer stabilen Molekel 
entspricht, durch Potenzentwicklung nach x gewonnen werden 


§ 2. Elektronenbewegung bei festen Kernen 


Setzt man in (12)x = 0, so bekommt man eine Differen- 
tialgleichung für die z, allein, in der die X, als Parameter 
vorkommen: 


1) Daß diese Auflösung im allgemeinen auf mehrdeutige Funk- 
tionen führen wird, spielt hier keine Rolle, ist aber physikalisch von 
Wichtigkeit; vgl. F. Hund, Ztschr. f. Phys. 43. S. 805. 1927. 
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(13) {H, (2,35; x) —Why=0. 


Sie stellt offenbar die Bewegung der Elektronen bei fest- 
gehaltenen Kernen dar. Wir nehmen an, daß dieses Eigen- 
wertproblem gelöst sei. Die Eigenwerte hängen nur von den 
Verbindungen £, der X, ab; denn man kann das Koordinaten- 
system mit dem System der Haupttrigheitsachsen zusammen- 
fallen lassen, wobei die X,=r,(£) werden. In diesem Achsen- 
system sind auch die Eigenfunktionen außer von den z, nur 
von den $&, abhängig; transformiert man aber auf beliebige 
raumfeste Achsen zurück, so treten daneben auch noch die 
9, auf. 

Wir bezeichnen den n-ten Eigenwert und die zugehörige 
normierte Eigenfunktion mit 


(14) W=V(8), vw=9 §, 9), 
so daB die Identitat 


(15) 89) — = 0 


gilt. Dabei nehmen wir an, daß /, ein einfacher Eigenwert 
ist. Tatsächlich wird dies nie der Fall sein, weil wegen der 
Gleichheit der Elektronen die von Heisenberg und Dirac 
entdeckte Resonanzentartung eintritt, bei zweiatomigen Molekeln 
überdies noch eine Entartung des Drehimpulses um die Achse. 
Aber da es sich uns hier nur um die Systematik des Näherungs- 
verfahrens handelt, wollen wir von diesen Entartungen absehen. 
Ihre Berücksichtigung würde auf Säkulargleichungen in den 
höheren Näherungen führen. 

Das wichtigste Ziel unserer Untersuchung ist der Nach- 
weis, daß die Funktion 7, (£) bei Bewegungen der Kerne die 
Rolle der potentiellen Energie spielt. Hierzu sind einige Hilfs- 
formeln nötig, die wir jetzt ableiten wollen. Es handelt sich 
darum, zu zeigen, _ die Matrix, die zu den Ableitungen des 


Operators H, 9) nach den £&, (also bei festgehaltenen 


2,35) gehört, ant ie Ableitungen der Funktion /, (§) zurück- 
geführt werden kann. 


Statt die Ableitungen nach den &, direkt zu bilden, er- 
setzen wir überall £, durch &, + x£, und differenzieren nach x; 
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der Koeffizient einer Potenz von x ist dann ein homogenes 
Polynom der ¢,, dessen Koeffizienten die Ableitungen nach 
den £, sind. So schreiben wir 


wo 

a) = 
b) 7% 


ö? Vi 
c) 


und in Weise 
HH, = Ayo + 2H? +...,; 
Man kann nun die Größen g,”, g,® nach pan Eigenfunk- 
tionen @,°(z; &,9) entwickeln; wir setzen: 


(17) 


(18) 


Dabei ist u’, ein as ae r-ten Grades in den 
&, z.B.) 
— 99% 


Die hier auftretenden nonin in denen dz das Volumen- 
element im Konfigurationsraum der Elektronen bedeutet, sind 
von der Orientierung des Kernsystems im Raume, also von 
den 9,, unabhängig; man kann sie etwa im System der Haupt- 
trägheitsachsen auswerten. 

Ist nun F irgendein auf die z, wirkender Operator, so 
wollen wir 


(21) [RP =F 
1) % ist die zu x konjugiert komplexe Zahl. 


(20) 
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das Matrixelement r-ter Ordnung von F nennen. Für r = 0 
geht es in das gewöhnliche Matrixelement 


(22) Poy = Fan = de 
über. Allgemein ist ae (19): 

Nun folgt aus (15) für x = 0 


Ferner ergeben sich durch Einsetzen der Reihen (16) und (18) 
in (15) zunächst die Identitäten: 


Multipliziert man diese mit %, und integriert über die x, so 
erhält man mit Rücksicht auf (24): 
a) u Ar + ar; =0, 


(26) b) (Fy — + (Hy — Vi aw 
+ (Han Onn! 0, 


Hieraus kann man der Reihe nach (H,')nw, (Haw, be 

/an 10S; /nn 

Von diesen Formeln werden wir nachher Gebrauch machen.’) 


§ 8. Aufstellung der Näherungsgleichungen 


Eine beliebige Konfiguration von Elektronen und Kernen 
läßt sich natürlich nicht mit einem allgemeinen Näherungs- 
verfahren behandeln. Wir wollen hier nur diejenigen Zustände 
ins Auge fassen, die einer stabilen Molekel entsprechen. Daher 
stellen wir folgende Frage an die Spitze: 


1) Das klassische Analogon zu der einfachsten Folgerung dieser 
Formeln, nämlich der aus (26a) für m =n’ folgenden Identität (H,"),„= Va” 
findet sich bei W. Pauli, Ann. d. Phys. 68. 8, 177. 1922; vgl. insbeson- 
dere 8 4, Formel (11). 


(28) 
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Gibt es ein Wertsystem der relativen Kernkoordinaten &, 
von der Art, daß die Eigenfunktionen w, des Energieoperators (6), 
soweit sie von den £, abhängen, nur in einer kleinen Nachbar- 
schaft dieses Wertsystems merklich von Null verschiedene 
Werte haben? 

Diese wellenmechanische Forderung entspricht offenbar der 
klassischen Bedingung, daß die Kerne nur kleine Schwingungen 
um eine Gleichgewichtslage ausführen; denn | , |? ist die Wahr- 
scheinlichkeit dafür, bei gegebener Energie eine bestimmte 
Konfiguration anzutreffen. 

) Als ungestörtes System betrachten wir also die Elektronen- 
bewegung für eine zunächst beliebig, aber fest gewählte Kern- 
konfiguration &,. Sodann entwickeln wir alle Größen nach 
kleinen Abweichungen von &,, die wir mit x¢, bezeichnen; wir 
setzen damit voraus, daß der „Schwingungsbereich“ mit x gegen 
Null geht, eine Annahme, die erst durch den Erfolg gerecht- 
fertigt werden kann. 

Sodann haben wir wie in § 2, (18) die Entwicklung 


4 a) H,°= H, (2,553 é), 
| | 


(28) 
c) = 05,08,’ 


und nach (11) wegen = 12. 
1 1 
(x, 3x)” Heo + Hoo) 


0 (1) (2) 
+ + Hye) + (Hy, + Hyg + 
0 ö° 
a) = Hy. (8 
anus 
7 


(29) 


wo 


(81) 


(32) 


(83) 


Die 
aus 


den 


(34) 


(35), 
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(9) He, = Hz, (8, 9, 


aH, 
(1) 


ar’ 


Die Argumente &; sind von jetzt an als Konstante zu betrachten. 


Der gesamte Energieoperator wird also: 
H= H,° + + + 
+ + + HE 
+ + Ho, + H + HS) +... 


folgenden Glieder haben alle dieselbe Form und entstehen 
dem Gliede mit «* der Reihe nach durch Erhöhung der 


oberen Indizes um 1. 


Wir entwickeln nun auch die gesuchte Eigenfunktion und 
Energieparameter nach x: 
y= 
+2: W®+.... 


Dann erhalten wir folgende Näherungsgleichungen: 


a) (H,° — =0, 

b) — = (W — HP) 

(HP — = — — 
+W® — Hw, 

d) (Hi? — = (WO — — — 
+(W® — H,® Hy) +(WV — y®, 

e) (H,° — W) = (WwW — — B® 
+(W® — — H°, — HD) yo 


+ (W® H,® + (W™ H,”) y®, 
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$ 4. Lösung der Näherungsgleichungen nullter und 
erster Ordnung: Kerngleichgewicht 
Die Gleichung nullter Näherung (35a) stellt die in 8 2 
diskutierte Elektronenbewegung bei festen Kernen dar. Aus 
der in § 2 eingeführten, zum Eigenwert 7° = V,, (&) gehörigen 
normierten Eigenlösung g,°(x; &, 9) erhält man die allgemeine 
Lösung in der Form 


(86) 
wo x,’ eine zunächst willkürliche Funktion der Argumente 
{, #, ist; diese muß man zur Verfügung haben, um die Lös- 


barkeit der folgenden Näherungsgleichungen zu erreichen. 

Die folgende Näherungsgleichung (35b) 
(37) (H,° W,) = “a H,”) y,° 
ist nämlich nur lösbar, wenn die rechte Seite auf w,° (bezüg- 
lich der Elektronenkoordinaten z,) orthogonal’) ist. Das gibt 
die Bedingung 
(38) = 
wo (H,®),„ das Diagonalelement des H,” bezüglich 
der x,, also nach (28b) eine homogene lineare Funktion der 
& ist. Diese soll nun nach (38) konstant sein, da 7,°(¢, 9) 
nicht identisch verschwinden kann, ohne daß dasselbe für w,° 
gilt. Daraus folgt 
(89) W%=0, (H,),, =0. 
Nach (26a) und (17) ist aber 


Also erhält man: 


(40) =. 


Die Fortsetzbarkeit unseres Näherungsverfahrens verlangt dem- 
nach, daß die relativen Kernkoordinaten &, nicht beliebig ge- 
wählt werden dürfen, sondern einem Extremwert der Elek- 
tronenenergie 7 (£) entsprechen müssen. Die Existenz eines 
solchen ist also die Bedingung für die Möglichkeit der Existenz 
der Molekel, ein Satz, der gewöhnlich als selbstverständlich 


1) Wir definieren die Orthogonalität zweier Funktionen f(x) und 
g(x) durch [f@)g (a) da = 0. 
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richtig angenommen wird. Wir werden nachher zeigen, daß & 2 
sich notwendig um ein Minimum handeln muß. Ker 
Die Funktion 7,°(¢, :#) bleibt zunächst unbestimmt. Setzt sich 
man in (37) W,°= =0 und w,° =7,°¢,! Ext 
ein, so erhält man für wy,‘ die Bestimmungsgleichung qua 
(41) (H,° Wa” H,® Fre 
Eine Lösung dieser ist nach (25a) gen 
gle 

YW, X n Pn „N 


wo 9,'” die durch (18) definierte Funktion (19a) ist. Die all- 
gemeine Lösung erhält man daraus durch Addition einer | "* 
Lösung g,° der homogenen Gleichung mit dem zunächst ur- all 
bestimmten Faktor 7,” (&, 9): 


(42) ap, ® Ku? D4 (47 
§ 5. Lösung der Näherungsgleichungen zweiter und De 
dritter Ordnung: Kernschwingungen zu: 


Wir gelangen nun zur Näherungsgleichung (35c), die nach J de 
Einsetzen der bereits gefundenen Lösungen der niederen § no 
Näherungen so lautet: 


Damit diese lösbar ist, muß wiederum die rechte Seite auf M 
g,,° orthogonal sein; das liefert mit den in § 2 eingeführten D 
Bezeichnungen wegen (39): in 
(Ho? + + = 0. 
Nun folgt aber aus (26b) für 0: 
Da ferner Hy. nach (30a) von den z, gar nicht abhängt, so 
- erhalten wir: ¥ 
(45) + = 0. ( 
Beachtet man hier die durch (17c), (30a) gegebene Bedeutung } 
von H,, und 7,®, so sieht man, daß (45) die Gleichung der 
harmonischen Kernschwingungen darstellt: ‘ 


0 6? 1 0°V, 
(46) (8 +3 why =o. 


BB 


en 
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Diese Gleichung zeigt, daß die Funktion /, (g) tatsächlich bis 
zu Gliedern 2. Ordnung die Rolle der potentiellen Energie der 
Kerne spielt. Für die Existenz einer stabilen Molekel ergibt 
sich nun die weitere Bedingung, daß der durch (40) bestimmte 
Extremalwert von 7, (£) ein Minimum sein muß; denn die 
quadratische Form 7, ® muß positiv definit sein, damit alle 
Freiheitsgrade ¢, stabil und Schwingungen um die Gleich- 
gewichtslage möglich sind. Bekanntlich ist die Schwingungs- 
gleichung (46) durch eine lineare Transformation der &, auf 
„Normalkoordinaten“ 1; separierbar. Ist o4,(£) die zum Eigen- 
wert W gehörige normierte Eigenlösung von (46), so ist die 
allgemeine Lösung 


Der Index soll also die Quantenzahlen der Kernschwingungen 
zusammenfassen. 9, (:#) ist eine zunächst unbestimmte Funktion 
der #,, deren Einführung für die Fortsetzbarkeit des Verfahrens 
notwendig ist. 

Bekanntlich ist 0), (£) eine lineare Kombination von Pro- 
dukten Hermitescher Orthogonalfunktionen der einzelnen 
Normalkoordinaten 7,; diese Funktionen haben die Eigenschaft, 
außerhalb des Gebietes der Librationsgrenzen der klassischen 
Mechanik sehr schnell (exponentiell) gegen Null abzufallen. 
Damit ist nachträglich unser Ansatz (£ +x£) gerechtfertigt, 
indem es sich zeigt, daß er tatsächlich auf Lösungen führt, 
die in bezug auf die &£-Oszillationen in Grenzen bleiben, die 
mit x verschwinden. Wir benutzen von den Hermiteschen 
Orthogonalfunktionen ferner die Eigenschaft, daß sie entweder 
gerade oder ungerade Funktionen ihres Arguments sind. 

Ist ® irgendein auf die ¢, wirkender Operator, so können 
wir die zugehörige Matrix 


(48) =[0 at 


ar) 


bilden, wo df das Volumenelement im Raume der £, ist. 


Um die Gleichung (43) zu lösen, setzen wir auf der rechten 
Seite nach (45) 
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ein; dann lautet (43) 
(Hr — V,°) pa” 

Ihre allgemeine Lösung ist 

(50) 

wo ee eine neue unbestimmte Funktion der ¢,, +, bedeutet; 


man überzeugt sich davon leicht auf Grund der Identitäten (25), 

Nun gehen wir zur Untersuchung der Näherungsgleichung 
3. Ordnung (35d) über, die nach Einsetzen der bereits bestimmten 
Größen so lautet: 


(H,° — V,°) = (W® — H® — He, — Hee) 
(51) + (#9 — — = (Lae Pr? + 


— Hy (xe, + +x 99%. 


Die rechte Seite können wir uns nach den g,° entwickelt 
denken; wir schreiben 


(49) 


wo 


nn’ nn’ Ans?’ 


dabei sind die F Operatoren nach £ und %, und zwar 
F®) — (H,®) , 
(54) | an 2 (2) (1) 
während wir von 7°) nur so viel auszusagen ur daß es 
eine homogene Funktion 3. Grades der ¢, und der 3% - ist. 
Damit (52) lösbar ist, muß 
sein; wegen (53) und (54a) bedeutet das 
(55) PE = (We — FO) , 
wo nach (54b) und (44) 


(58 


ist. ( 
inhot 
die 1 
Integ 
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Abe: 
ist, 
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|: 
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wo 
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ist, (55) ist also die zur Schwingungsgleichung (45) gehörende 
inhomogene Gleichung; da (45) die zum Eigenwert W ge- 
hörige normierte Lösung o° hat, so ist (55) nur lösbar, wenn 
die mit a multiplizierte rechte Seite ein verschwindendes 
Integral über den ¢-Raum hat. Das gibt wegen (47b) eine 
Differentialgleichung für 9° (9): 
Pe) _ = 0. 
Aber da, wie wir sahen, 7° ungerade in den £, und =: 


ist, so muß das Diagonalelement der zugehörigen ¢-Matrix 
verschwinden. Denn transformiert man auf die Normalkoordi- 
naten 7,, so wird o° eine Summe von Produkten Hermite- 


scher Orthogonalfunktionen, F‘ ein Polynom ungerader Ord- 
nung in den n, und oa so daB jedes Glied mindestens eines 
der 7, oder 7 in ungerader Potenz enthält; daher ver- 


schwindet bei der ¢-Matrixbildung jedes Glied einzeln. Es 
folgt also 


(56) W®=0, 
und 0, bleibt vorlaufig weiter unbestimmt. 
Nunmehr läßt sich (55) lösen: 


wo gs" folgender Operator bezüglich der +, ist: 
w® _ w® 


Endlich wird die Lösung von (52) 


und das hat nach (53) die Bi 
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wo 


(3,2) 

(61) _ 
nn v? v° 
n n 


gesetzt ist. Nach (54) ist G° eine Zahl, G®”) ein Differential 
operator bezüglich der ¢,, G*” ein Operator begüglich der { 
und 
Nach § 2, (26a) ist 

ay 

a) 


= Fn 
also wird 


3 1) , (2) (3,2) „(U 0 (3,3) „0 0 
(62) Ane + > Pr + Gin’ 


§ 6. Lösung der Näherungsgleichungen 
vierter und höherer Ordnung: Rotationen und Koppelungseffekte 


Die Näherungsgleichungen 4. Ordnung (35e) lauten nach 
Einsetzen der bereits bestimmten Größen: 


— = 
= (9% — — — 2, - Be) 

(1) (0) 0 9 (1) 
(63) = Hy, + H;,) (Zu + ) 


(2) (2) 0 (2) 0 (1) (1) 0 _® 
Hi, + H,.) (xn + 2,9, + 


nn 


Wir entwickeln wieder die rechte Seite nach den g,°: 
= — D 


nn’ 
wo 


(4) 4,2) (2) 4,3) (1) 4,4) 0. 
hier ist 

4,2) 0 (2) , (2 (1)\(1) 
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und das ist mit 7°” (54b) identisch, während 7“ ungerader 


nn 
Ordnung in den ¢,, & ist, PC) aber gerader Ordnung. Die 
Integrabilitätsbedingung von (64) lautet: 
0 4) 
das bedeutet nach (65) 
nn ns nn n nn ns 
Die linke Seite stimmt wegen (66) wieder mit der der 
Schwingungsgleichung (45) überein. Die rechte Seite muß also 
auf o? normal sein. Setzt man die Werte von 7°) und 7°) 
aus (47b) und (57) ein und benutzt das Symbol 


68) 


so erhält man 


Diese Gleichung bestimmt endlich die Funktion 0), ($), also 
die Bewegung des Systems der Hauptträgheitsachsen im 
Raume, die Translationen und Rotationen. Das Hauptglied 
des Operators in (69) ist dasjenige, welches die zweiten Ab- 
leitungen nach den », enthält; ein Blick auf (63) zeigt, daß 
es aus HY, entsteht, ihm entspricht also in das 
Glied 

(10) 

wo an die Stelle der Punkte die Funktion einzusetzen ist, auf 
die der Operator wirkt. Da der Operator (70) von den ¢, un- 
abhängig ist, so sind die Diagonalglieder der zugehörigen 
s-Matrix mit ihm selbst identisch. Physikalisch bedeutet die 
Tatsache, daß statt des einfachen Operators H, » der kompli- 
ziertere (H},), auftritt, eine Koppelung der Kreiselbewegung 
der Kerne mit der Elektronenbewegung. Es handelt sich, wie 


wir nachher am Beispiel der zweiatomigen Molekel sehen 
Annalen der Physik. IV. Folge. 84. 31 
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werden, um diejenigen Effekte, die Kramers und Pauli)) 
durch die Annahme eines im Kreisel eingebauten „Schwung- 
rades“ darzustellen versucht haben. Sodann gibt es in (69) 
Terme, die sich aus dem Operator H,, ableiten; sie ent- 
sprechen einer Koppelung der Kreiselbewegung mit Dreh- 
impulsen, die infolge der Kernschwingungen entstehen. Endlich 
sind Glieder da, die auf die %, nicht wirken; sie geben einen 
Zusatz der Ordnung x* zur Schwingungsenergie. 

Da die Translationen sich immer in trivialer Weise ab- 
separieren lassen, betrachten wir nun die Rotationen. Ist r 
die zugehörige Quantenzahl, so hat man für die Lösung von (70) 


(71) we yw. 0, (9) 


ner? 
zu schreiben. Sodann kann man (67) lösen und schließlich 
auch (64). Doch hat es keinen Sinn, die Formeln explizite 
anzuschreiben. 

Es ist klar, daß das Verfahren sich fortsetzen läßt, ohne 
das sich etwas prinzipiell Neues dabei ergibt. Die höheren 
Näherungen beschreiben die Koppelungen zwischen Rotationen, 
Schwingungen und Elektronenbewegungen. Neue Quantenzahlen 
außer den bereits eingeführten treten nicht mehr auf. 

Wir fassen nun die charakteristischen Züge der gefundenen 
Lösung zusammen. Das augenfälligste Ergebnis ist, daß zur 
vollständigen Bestimmung der Eigenfunktionen auch nur in 
nullter Ordnung die Lösung der Näherungs-Differential- 
gleichungen bis zur 4. Ordnung notwendig ist; man hat nämlich 


72) yy, 


wo g,° die Eigenfunktion der Elektronenbewegung bei fest- 
gehaltenen Kernen, a , die der Kernschwingungen und 0) a 
die der Rotationen sind. Dabei sind die Schwingungskoordı- 
naten ¢; von einer Gleichgewichtslage &, aus zu zählen, die 
durch die Forderung bestimmt ist, daß in ihr die Energie der 
Elektronenbewegung 7, (£) ein Minimum ist. Die Bestimmung 
der drei Funktionen p,°, o° , 0° 


ns? nay iefert zugleich die Energie 
bis zur 4. Ordnung einschlietlich: 


1) H. A. Kramers und W. Pauli, Ztschr. f. Phys. 13. S. 343, 
351. 1923. 
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dabei ist V,° der Minimalwert der Elektronenenergie, der die 
betrachtete ruhende Molekel kennzeichnet, W die Energie 
der Kernschwingungen, W) enthält (neben Zusatzgliedern zur 
Schwingungsenergie) die Energie der Rotationen. In dieser 
Nähernng (bis x*) erscheinen also die drei Grundtypen der 
Bewegung „separiert“; will man ihre Koppelungen untersuchen, 
muß man höhere Potenzen von x berücksichtigen. 5 

Auf Grund von (72) lassen sich Aussagen über die Über- 
gangswahrscheinlichkeiten (Intensitäten von Bandenlinien) 
machen. 

Das elektrische Moment der Molekel M setzt sich aus 
einem Anteil der Kerne ® und einem Anteil der Elektronen p 
zusammen; eine Komponente sei 

{ %, > e, x, 2 
(74) M, = 
k 
Hieraus bilde man der Reihe nach die Matrixelemente bezüg- 
lich der z,, &, 9,; zunächst ist 
(75) (Pen, = {Pe g, az 
eine Funktion der ¢,, o;, sodann sind 


(8) = [Pan ls (Bons = 
nur noch Funktionen der a, endlich 

| (Poin = f Pons On a'r dt, 

\ Bln or = [Bons 

numerische Konstanten, die die Ausstrahlung und die Uber- 


gangswahrscheinlichkeit für nsr—>n's’r' bestimmen. Diese 
schrittweise Bildung kann man so interpretieren: Zu jedem Elek- 
tronensprung n—>n' gehört ein virtueller Oszillator mit dem 


Moment (p,),; aus diesem gewinnt man die Matrix (p,), 4? die 
n’ 


m) 


einem System von Schwingungsbanden (Übergänge s—> s’) ent- 

spricht, nach einer von der üblichen etwas abweichenden Regel, 

indem man eine Eigenfunktion des untern, eine des obern 

Elektronenniveaus benutzt [Formel (76). Das Entsprechende 
s1* 
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wiederholt sich bei den Linien einer Bande (Übergänge r —> r', 
Darin sind die zuerst von Franck!) angegebenen, von Condon* 
weiter ausgeführten Abschätzungen der Intensität von Schwin- 
gungsbanden enthalten. Diese wird nämlich bestimmt durch 
den Verlauf der Funktionen 7 (&), 7„(&); nur in der Um 
gebung ihrer Minima sind die zugehörigen Eigenfunktionen ot | 
co}, , merklich von Null verschieden, ihr Produkt also nur dann, 
wenn diese Bereiche übereinander greifen. Wenn die Funk- 
tion 7, (&) bei einem Elektronensprung n—>n’ nur wenig in 
ihrem Verlauf verändert wird, werden die Banden mit kleiner 
Änderung von s intensiv sein; sobald aber /, (&) sich bei dem 
Sprunge n—> n’ stark wong, wird ein Übereinandergreifen 
der Intervalle, in denen of, und o},, nicht verschwinden, nur 
bei größeren Differenzen s — s’ möglich sein. Diese Verhält 
nisse sind von Condon quantitativ diskutiert worden. Ähnliche 
Betrachtungen müssen sich mutatis mutandis auch für die 
Rotationen durchführen lassen. 


§ 7%. Spezialfall der zweiatomigen Molekeln 


Als Beispiel wollen wir die zweiatomigen Molekeln kurz 
behandeln. Bei diesen hat man außer der Resonanzentartung, 
die von der Gleichheit der Elektronen herrührt, noch eine 
weitere Entartung, indem zu jedem Energiewert zwei Be 
wegungen gehören, bei denen der Drehimpuls um die Kern- 
verbindungslinie entgegengesetzt gleich ist. Da es sich uns 
hier aber nicht um die wirkliche Diskussion der Feinstruktur 
der Banden handelt, wollen wir auch von dieser Entartung 
absehen; wir beschränken also unsere Betrachtung auf Fälle, 
wo der Drehimpuls um die Achse verschwindet bzw. die 
Energie der Elektronenbewegung nicht oder nur wenig von 
ihm abhängt. 

Bei zwei Kernen haben wir nur eine &-Koordinate, on 
Kernabstand £, und 5 #-Koordinaten, nämlich die Koordinaten 
des Schwerpunktes X,, F,, Z, und die Polarkoordinaten 9, o 
der Kernverbindungslinie (Achse). 


1) J. Franck, Trans. Faraday Soc. (1925). 


2) E. Conden, Phys. Rev. 28. $. 1182. 1926; Proc. Nat. Acad. 13 
8. 462. 1927. 
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Die kinetische Be Kerne wird 
+ 


wo 
(M, + M,)® 
(78) 
und 
6° ö 6° 
[4 = + ay, 
(80) | 1 ö 
I, = sin?’ OF (sin 38) 
gesetzt ist. 
Wir haben also: 
h? 0° 
H,=- tam 
ö 
(81) H,, + 
Ersetzt man hier & durch § +f und ami nach x, so 
erhält man: 
(82) 
HY? =0, p=1,2, 
2u 
(83 
) | Sa? m 6c 
h? 
Hy, = Er (4 5 4,) 
84 (1) 2u 
(84) saw 


(85) 
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Der Kernabstand £ bestimmt sich durch die Gleichung 

= 
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Die Gleichung der Kernschwingungen (45) lautet hier: 


[St 


Setzt man nun 


8n? m (2) _ In? ” 
4 
n= Yb, 


so hat man’) 
Die Eigenwerte sind 
ye (s = 0, 1, 2,..,), 
die Eigenfunktionen 


wo H, das s-te Hermitesche Polynom ist. 
Die Schwingungsenergie ist also: 


hi 


m 
=(s+ 7, 
oder 
(88) = (s+ 
dabei ist 
1 1 1 1 
(89) = (+ az) = 


die Frequenz des Oszillators. 

Wir stellen nun die Gleichung (69) für die Rotationen 
auf, verzichten aber auf ausführliche Angabe der Korrektur 
an der Schwingungsenergie. Da hier H, ze nach (81) die Ab- 
Ieitungen nach den 9, gar nicht enthält, so kommt für uns 

in (69) nur der Term Uz (Ae 5)n in Betracht; alle übrigen Glieder 
fassen wir in die Konstante C,, zusammen. Dann lautet die 
Rotationsgleichung (69): 


(90) + Ons — Wl 0°, = 0. 


1) 8.E. Schrödinger, Ann. d. Phys. 79. $. 361, $ 3. 1926. 
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Indem wir in H} » den translatorischen Anteil weglassen, haben 
wir nach (70) und (84) für - beliebige Funktion (9): 


As (py f)dz 


und nach (80) 

= 9, 4of +1409, +2 

Daher wird 

f = - nd 


ı af dg? 
+2 Go + ar). 


Indem wir 4, in - ae 


1 0? 
A, = t+ BOSS + 


schreiben, sehen wir, daB es angebracht ist, folgende Ab- 
kürzungen einzuführen: 


Diese Größen sind die der von P,, 
P,, bzw. p,”, p,,” (von einem Fakto 
bedeuten also die Mittelwerte der zu i ‘welled Winkel 
gehörigen Drehimpulse der Elektronenbewegung, bzw. die Mittel- 


werte ihrer Quadrate. Dann lautet die Gleichung (90) aus- 
führlich: 


(91) 


Kor + 26. 39 +) + 989 + 6,) 


02) + (Gor + + OD) 


— = 0 
Sie ist der von ae und Pauli aufgestellten Gleichung 
für den Rotator mit eingebautem Schwungrad sehr ähnlich; 
der Unterschied besteht im wesentlichen darin, daß bei diesen 
Autoren statt der, Quadratmittel ©,?, 2,? die Quadrate der 


Mittelwerte, auftreten. 
Übrigens läßt sich die Abhängigkeit der Größen (91) von 
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den Winkeln 3, wm durch eine elementare Betrachtung fest- 
stellen, wenn man annimmt daß für diesen Zweck die Diagonal- 
elemente der quantenmechanischen Matrizen durch die ent- 
sprechenden klassischen Mittelwerte ersetzt werden dürfe. Die 
Bewegung des Drehimpulsvektors der Elektronenbewegung kann 
man zerlegen in eine unregelmäßige Schwankung ohne mittlere 
Rotation und eine überlagerte gleichférmige Rotation um die 
Molekelachse. Die Schwankung ersetzen wir im Mittel durch 
einen festen Vektor; dieser rotiert dann gleichförmig um die 
Achse. Man hat also dasselbe Verhalten wie bei einem sym- 
metrischen Kreisel, dessen Drehimpuls bezogen auf ein im 
Kreisel festes Bezugssystem die Komponenten Z, M, N haben 
möge. Aus diesen lassen sich die Komponenten des Dreh- 
impulses, die zu den Winkeln 9, w gehören, so ausdrücken‘): 
0 = Lcosy — Msiny, 
2= Lsindsiny + Msind cosy + Ncos#; 

dabei ist y der Winkel der Eigenrotation um die Achse. Durch 
Mitteilung über 7 folgt daraus: 

6 =0, 2 = Noose, 

=4(L?+ M%), 2? =4(L? + + N2cos?#. 
Hier hat man N mit der Quantenzahl p zu identifizieren, die 
den Drehimpuls um die Achse gibt, und 4(Z? + 72), 4.N? mit 
den Quadratmitteln p*, p> des gesamten Elektronendreh- 
impulses senkrecht und parallel zur Achse; da N konstant ist, 
gilt p> =p’. Man hat also schließlich 


6, =0, 2, = pcos#, 
(98) =p?, =p? + p? cos*#. 
Dieses Resultat bedarf aber natiirlich einer strengen quanten- 


mechanischen Nachprüfung; vermutlich ist p? durch p(p + 1) 
zu ersetzen. 


Durch das Eigenwertproblem (92) wird die Größe 
an W™ gleich einer numerischen Funktion der Rotations- 
quantenzahl r, etwa gleich g, ‚(r); für die Rotationsenergie folgt also 


ner Ins 


1) Vgl. etwa: F. Klein und A. Sommerfeld, Theorie des Kreisels 
1. S. 108. 
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das Trägheitsmoment der Kerne im Gleichgewichte ist. 

Eine Diskussion der höheren Näherungen hat ohne Ein- 
gehen auf die Entartungen keinen Sinn und soll hier unter- 
bleiben. 

Dagegen wollen wir noch kurz zeigen, daß man die zwei- 
atomige Molekel durch ein ganz anderes Störungsverfahren 
behandeln kann, dessen klassisches Analogon von Born und 
Hückel!) durchgeführt worden ist. Dabei wird als ungestörte 
Bewegung nicht das Elektronensystem bei ruhenden Kernen 
genommen, sondern bei gleichförmiger Rotation der Kerne. 


$ 8. Unabhängige Behandlung der zweiatomigen Molekeln 
Wir gehen auf die Gleichung (12) zurück und schreiben 
sie unter Einsetzen von (11): 
{H, + x*(H, st Ho + Hos) - W\y=0. 
Bei den anieeeiiiinn “Molekeln besteht nun die Besonderheit, 
daß H,, von & überhaupt nicht abhängt. In diesem Falle 
führt die Methode zum Ziele, Translationen und Rotationen 


sogleich abzuseparieren. Man hat nämlich nach (81) (unter 
Weglassung = 


(96) wat #40) - W}y=0. 
Setzt man nun 
(97) = 8); 
wo Y, eine Kugelfunktion r-ter Ordnung ist, also der Gleichung 

genügt, so behält man für #, die Bedingung 

1 

(98) Im — (Sp - 
Nun ersetzen wir wieder durch &+xZ, fragen also nach 
Schwingungen um den Zustand gleichförmiger Rotation. Be- 
zeichnen wir die Energie dieses Zustandes mit 

2 2 2 


1) M. Born und E. Hückel, Phys. Ztschr. 24. S. 1. 1923. 
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und setzen 

(100) W=E+R, 

so wird die Gleichung (98) 
(101) (H°+ +...—-E)¥ =0, 
wo 


H’ = H,°, 
HY — HY +R, 
1 


H® = H,® + OR” — eu 2 


Of’ 


H,°®, H,®, .... bedeuten die früher ebenso bezeichneten 
Operatoren. Es gelten alle Formeln von $ 2 unverändert. 
Die Näherungsgleichungen lauten 
a) =0, 
(108) b) (H® yp (Eo H‘?) po, 
c) (H° E°) (E H®) yo + (E® H®) pn, 

Die erste hat die Lösung 
wo V,(&), »,°(x; §) die früher eingeführten Funktionen und 
s,,(£) vorläufig willkürlich ist. Die Integrabilitätsbedingung 
für (103b) lautet 

(ET - Q=0. 
Nun ist nach $ 2, (26a) 

Ho) + ER = = CA, 

Daher folgt, wie früher (§ 4), daß 
(105) E™=0, =0. 


Diese Bedingung besagt offenbar, daB bei der ungestérten 
Rotation Gleichgewicht herrschen muß zwischen der Zentrifugal- 
kraft und der quasielastischen Kraft, die infolge der Elektronen- 
bewegung einer Verrückung der Kerne widersteht. Die Zentri- 
fugalkraft ist nämlich 
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1 1\ A r(r +1) 
(Ge + aig) (Ge + ag) 
wo der quantenmechanische Wert + Vr(r + 1) für den Dreh- 


impuls p, eingesetzt ist, und das stimmt nach (99) und (95) 
mit 2’ überein. 


Aus der Beziehung (105) hat man den Gleichgewichts- 
abstand € zu berechnen, der noch von der Rotationsquanten- 
zahl r abhängt. Für kleine Werte der Rotationsenergie 2 kann 
man £ nach Potenzen von 

2 
entwickeln; man erhält!) 

1 3 A 
00) 


Da 8 von der Ordnung x‘ ist, wird man bei systematischem 
Vorgehen nur so viele Glieder dieser Reihe gebrauchen, als 
dem Grade der Approximation beim Störungsverfahren ent- 
spricht. 

Indem wir uns nun diesem wieder zuwenden, wollen wir 
uns kurz fassen, da die Methode dieselbe ist wie früher, nur 
vereinfacht durch die Vorwegnahme der Rotationen. Die Lösung 
von (103b) lautet igang (42) 


und die Integrabilitätsbedingung. von (102) liefert 
Das ist aber die RENNER 
» (2) 


Es wird also, wie in § 7, 


(110) = (5+ 
wo die Frequenz 


1) Man kann diese Formel einfach aus der zitierten Arbeit von 


Born und Hückel entnehmen. 
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noch von der Rotationsquantenzahl r abhängt, die in A vor. 
kommt. 
Ferner wird, wie in § 7, 


(112) o =e 
mit 

4 2 
(112a) n=CYb, b= +R’). 


Das Verfahren läßt sich in der bekannten Weise fortsetzen, 
Man findet E® = 0, während E® außer der Abweichung vom 
harmonischen Schwingungsgesetz auch eine Koppelung mit der 
Elektronenbewegung enthält. Doch würde die ausführliche 
Angabe der Formeln aus dem Rahmen dieser Abhandlung 
fallen, die nur das Prinzip der Entwicklung zur Darstellung 
bringen soll; auch hat die Durchrechnung der höheren Nähe- 
rungen nur Sinn, wenn gleichzeitig die Entartungen berück- 
sichtigt werden. 


H, (4) 


(Eingegangen 25. August 1927) 
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2. Wirbelströme im elliptischen Zylinder; 
von M, J. O. Strutt 


Wir denken uns einen sehr langen zylindrischen Leiter 
yon elliptischem Querschnitt in einem longitudinalen elek- 
trischen, bzw. in einem longitudinalen magnetischen Wechsel- 
feld. Diese beiden Fälle entsprechen der „Stromverdrängung“, 
bzw. der „magnetischen Feldverdrängung.“ Der Einfachheit 
halber nehmen wir das äußere Wechselfeld jedesmal homogen 
an. Physikalisch läßt sich dieser Zustand, wenn die Wellen- 
länge im Äther der elektromagnetischen Schwingung groß 
gegenüber den Querschnittsabmessungen des Zylinders ist, für 
das longitudinale elektrische Wechselfeld näherungsweise reali- 
sieren, indem man den Leiter zwischen zwei leitenden Platten 
ausspannt, auf die eine Wechselspannung wirkt. Für das 
longitudinale magnetische Wechselfeld kann man den Leiter 
in einer langen Magnetisierungsspule gelegen denken, durch 
die ein Wechselstrom fließt. 

Die Lösung des Problems ist in beiden Fällen formal 
gleich. Nur wenn man zu den numerischen Ergebnissen vor- 
dringt, ergeben sich Unterschiede, da man im Falle des longi- 
tudinalen elektrischen Feldes das Verhältnis der Impedanz 
zum Gleichstromwiderstand, im Falle des longitudinalen magne- 
tischen Feldes aber die im Leiter enthaltene Joulesche und 
magnetische Energie zu kennen wünscht. 

Wir werden dementsprechend bei der Lösung als ab- 
hängige Veränderliche eine Größe u einführen, die in den 
beiden Fällen bzw. mit der elektrischen Leitungsstromdichte 8 
und mit der magnetischen Feldstärke H zu identifizieren ist. 


1. Das Randwertproblem 


Die Leiterachse fällt mit der z-Achse zusammen. Wir 
setzen 


ler 
he 
Dg 
ng 
k- 
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Dann genügt u im Leiter der Differentialgleichung: 
um 
v = spez. Widerstand, 
u = Permeabilität, 
w = Kreisfrequenz, 
i=y-1l, 
V, = Lichtgeschwindigkeit. 

Es sind die Gaussschen Einheiten benutzt und es ist der 
Verschiebungsstrom gegenüber dem Leitungsstrom vernach- 
lässigt. 

Die Lage der Achsen und die Abmessungen der Quer- 
schnittsellipse sind der Fig. 1 zu entnehmen. 


Führen wir elliptische Koordinaten ein vermittelst der 
Gleichungen 
(2) { zr=c-cosy Gof &, 
y=c-sing Sing, 
(c = lineare Exzentrizität), 
wobei die Koordinaten gebildet werden von den Scharen kon- 
fokaler Ellipsen, bzw. Hyperbeln: 


+ 

(3) : 
cos? sin®n 


(4) 
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so schreibt sich die Gleichung (1): 


0° 0°? 
ger u (Gof 2g — cos 2] 
mit 
2h = 


Zur Lésung setzen wir 
u = V(§)- W(n), 


wobei V und W bzw. den Differentialgleichungen: 


(5) + W+(R + 2h? cos 2) = 0 
und 

6) V-{R + 2h? Coj 2H =0 
geniigen. 


Hierbei ist # eine zunächst ganz willkürliche Konstante. 
Da (6) aus (5) hervorgeht, wenn man 7 durch 7 £ ersetzt, werden 
wir uns vorerst nur mit der Lösung von (5) beschäftigen. 


Die Randbedingung ist dadurch gegeben, daß für 


&=& (Querschnittsumfang) 
gelten muß: 
u = Konstante. 


2. Die Lösungen der Mathieuschen Differentialgleichung. 


Die Differentialgleichung(5) wurde zuerst von E. Mathieu!) 
betrachtet im Zusammenhang mit den Schwingungen einer 
elliptischen Membran. 

Es kann sich in diesem Abschnitt nur darum handeln, 
kurz die wichtigsten Eigenschaften der Gleichung (5)?) zu über- 
blicken. Mathematisch werden wir dabei kaum etwas wesent- 
lich Neues vorzubringen haben. 


1) E. Mathieu, Cours de Mathématique physique S. 122. 1873. 

2) E.T. Whittaker u. G.N. Watson, Modern Analysis S. 404. 
1920. P. Humbert, Fonctions de Lamé et de Mathieu 1926, Gauthier- 
Villars. F. Tisserand, Traité de Mécanique Céleste III. S. 1. 1894. 
H. Poincaré, Méthodes nouvelles de la Mécanique Céleste II. S. 228. 
Hm. B.v.d. Pol jr. bin ich für viele Literaturangaben zu großem Dank 
verpflichtet. 
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Es sei f(n) eine Lösung von (5), so ist 
f(a +2) 
ebenfalls eine Lösung und es können zwei Fälle eintreten. 
a) fm 
Wir behaupten, daß im Falle F(n) = Lösung von (5), gilt 


v = Konstante. 
Im Falle b ist diese Behauptung bereits von f selber 
erfüllt. 
Im Falle a) gelte: 
Dann behauptet unser Satz (v, = Konstante): 


fa + 2a) fm 


mit 


oder: 


A-f(n)+ 
folglich 


1 _f@+F(-n 


v Fo) 
Setzen wir ' 
Fin) = a. z{n) + -y(n), 
wobei 
xz = gerade Funktion, 
y = ungerade Funktion, 
so gilt 


F(x) = a-2(n) + b-y(z), 
F(— 2) = a-2(x) — b-y(m), 
2 F(o) =a, 
folglich 
y+ = 2-r(n). 
Setzen wir weiter 
» = erp(ika), 


so ist 


coska = x(n). 


4 
die 
eir 
ist 
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Be 
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Es können nun drei Fälle eintreten: 


a) |coska|<1, 
b) |coska| =1, 
|coska|>1. 


Im Falle a) sind die Lösungen von „stabilem Typus“, 
d. h. die Ausschläge sind stets endlich. 

Im Falle c) wachsen die Ausschläge über alle Grenzen, 
die Lösungen sind „labil“. 


Im Falle b) sind die Lösungen ebenfalls labil.') Da x(z) 
eine Funktion der Parameter & und A? in der Gleichung (5) 
ist, entspricht jedem Wert von k, oder besser von coskn eine 
bestimmte Beziehung zwischen R und h?. Wir können diese 
Beziehungen als Linien darstellen in einer Ebene, mit 2h? 
als Ordinaten und R als Abszissen. Die Kurven überdecken 
diese Ebene überall dicht. Einige unter ihnen entsprechen 
labilen, andere stabilen Lösungen. Die Kurven, welche 
hervorgehen aus 

|coska| = |2(m)|=1, 


bilden die Grenzen zwischen den labilen und stabilen Lösungen. 
Einige dieser Kurven sind in der Figur 2 gezeichnet worden 
nach Reihen von Mathieu.) Man kann aus diesen Reihen 
leicht die Beziehungen nachweisen, die Poincaré’) in einer 
schönen Untersuchung über die Berührung dieser Kurven für 
2h? = 0 gefunden hat auf einem Wege, der sich an frühere 
Berechnungen von Tisserand‘) anschließt. Die periodischen 
Lösungen, welche der zuletztgenannten Beziehung entsprechen, 
haben die Periode 27. 

Bei der Berechnung von Lösungen der Gleichung (5) geht 
man aus vom Satze, daß sich jede Lösung nach Potenzen 
von h? entwickeln lassen muß. 


1) An anderer Stelle kommen wir auf die Bedeutung dieser Tat- 
sachen für gewisse mechanische und elektrische Probleme zurück. 
2) a.a.0. 8.132, u. w. 
8) a.2.0. S. 243. 
4) a.2.0. 8.8. 
Annalen der Physik, IV, Folge, 84, 
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Wir nennen die Lösungen, welche sich für 
h?=0 
reduzieren auf . 
cosgyn und sin gy, 
R (h?),->0 = 7 ’ 
Ce,(h’, n) und Se,(h?, m), 
und bezeichnen sie als Mathieusche Funktionen im allgemeinen 
Sinn. 
75 T 
| 


LABIL 
h? 
40 


| LABIL 


wobei gilt: 


bzw. 


STABIL STABIL 


- 
05 R 05 10 


Fig. 2 
Whittaker’) betrachtet lediglich solche Lösungen, für die 
g = ganze Zahl, 
und folglich 
|coska| = 1, 
und bezeichnet diese mit dem Namen Mathieusche Funk- 
tionen. Wir werden die Funktionen bezeichnen als Mathieusche 
Funktionen im engeren Sinn. 
Allgemeine Funktionen von Mathieu betrachtete bereits 
E. G.C. Poole.*) Für sie gilt 
|coska| <1 


1) a.a. 0. S. 410. 
2) Proc. Lond. math. Soc. 20. S. 374. 1922. 
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Sämtliche Mathieuschen Funktionen liegen im stabilen 
Gebiete, bzw. auf der Grenze zwischen den stabilen und labilen 
Gebieten der (2, h?)-Ebene. 

Für unseren Zweck ist es erforderlich, daß die lösenden 
Funktionen eindeutig vom Ort abhängen. Es kommen also 
nur die Mathieuschen Funktionen im engeren Sinne in Betracht. 


3. Mathieus Funktionen im engeren Sinne 
Reihen für sämtliche Mathieuschen Funktionen findet 
man bei Mathieu!) und Whittaker.?) 
Im stabilen Gebiete, d.h. für 


|coska| <1, 
geht von jedem Punkte der R-Achse nur eine Kurve 


f(R,A?) = 0 
aus, für die coska einen konstanten Wert besitzt und für 
jeden Punkt dieser Kurve gibt es zugleich eine Lösung Ce, (h?, 1) 
und eine Lösung Se, (h,,n), mit un Index g. 

Von den Punkten 1, 4, 9... der R-Achse gehen jedoch 

je 2 Kurven 

h?) =0 
aus, die der Beziehung 

|coska| = 1 


entsprechen. Längs der links gelegenen gibt es eine Lösung Se, 
md längs der rechts gelegenen eine Lösung Ce, Da nun 
aber das allgemeine Integral in jedem Punkte der (R, h?)-Ebene 
sich aus 2 Partikularintegralen zusammensetzt, gibt es längs 
jeder dieser Kurven noch nicht-periodische Lösungen®), die 
aber für uns, da sie keine eindeutige Funktionen des Ortes 
in unserer Ellipse sind, nicht in Betracht kommen. 

Die Mathieuschen Funktionen im engeren Sinne genügen 
den Orthogonalitätsrelationen 9): 


M + m -dn=0 (m + n) 


1) E. Mathieu, a. a. O. S. 130. 

2) E. T. Whittaker, a. a. O. S. 411. 

3) B. Sieger, Diss. Wiirzburg, S. 19. 1908. 
4) E.T. Whittaker, a. a. O. 8. 411. 
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(8) Se, 2) Se, = 0, (m +) 


(9) (h?, Ce, (h*,4)-dn 


In der Form, wie ihre Reihenentwicklungen vorliegen, 
sind Mathieus Funktionen nicht normiert, d.h. es gilt: 


m)? dy = Konstante, 


= 


f - dy = Konstante. 
Weiter sind sämtliche Ableitungen der Mathieuschen 
Funktionen, gemittelt über eine Periode, gleich Null: 


=0 p>o 


wobei p die p-fache Differentiation nach 4 bedeutet. In dem 
Intervall -=n=n gemittelt, sind alle Funktionen Ce w- 
gerader Ordnung und sämtliche Funktionen Se gleich Null: 


- dn=0, 


Jf Sen dn = 0. 


Lösungen der Gleichung (6) erhalten wir, indem wir 7 durch 
i£ ersetzen, aus denjenigen der Gleichung (5). Wir bezeichnen 
sie mit 

und Se, (h?, &). 

Bei gleichem Index g gehört zu jeder dieser Funktionen 

derselbe Wert 

R=fW), 
wie zu den entsprechenden Funktionen 
Ce,(h’,n) und Se,(h?, n). 


die 


(10) 


| 


Diffe 
habe 
und 
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4. Die Lösung des Randwertproblems 


Nach dem Vorhergehenden können wir eine Lösung der 
Differentialgleichung (4) ansetzen in der Form: 


> A, + en &) Ce„(h}, 2) 


+ Sen (h, &) 7). 


Zur Bestimmung der Konstanten 4, und B, haben wir 
die Bedingung, daß u am Leiterumfang (£ = &,) den Wert X 
haben muß: 


ent, = K= Cen 


+ SB, &,) (42, 7). 


Die Konstanten A, und B, erhalten wir, wenn wir links 
und rechts des Gleichheitszeichens bzw. mit Ce, (h?,) und mit 
Se,(h?,7) multiplizieren und sodann von —z bis m mitteln. 
Wegen der im vorigen Abschnitt erwähnten Beziehungen fallen 
alle B, und die A, mit ungeradem Index heraus und es ent- 
stehen: 


Fi C (h?, dn 


An, = K- n=0,1,2... 


STC ean (lt, -dn 
Die vollständige Lösung unseres Randwertproblems lautet 
somit: 
ay 
u=K. > 
(10) 0 
Ce,» (h?, &,) + dy 


Ce, „(h?, &) . Ce,, (h?, 1)» 
mit n= 0, 1, 2, 8... 


5. Reihenentwicklungen fiir die benutzten Funktionen 


Die von uns benutzten Funktionen Ce,, (h?, 7) sind iden- 
tisch mit denjenigen, die Heine!) als Funktionen der ersten 
Klasse bezeichnet. 


1) E. Heine, Handbuch der Kugelfunktionen, I. S. 401. 


| 
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Ihre Reihenentwicklung lautet: 


(11) (h®, n) = +d a,- Cos2nn, 
1 


wobei die Koeffizienten « Funktionen von Ah? sind. 

Bei der Berechnung der Reihen geht man so vor, daß 
man nach Potenzen von Ah? entwickelt. Bis zur 10. Potenz 
von A sind diese Reihen bekannt. 

Sie lauten): 


4 
Ce, (0, n) = cos2y +  cosdn + + 00867 
cos 87 43 cos 47 5 


cos 107 287-cos 67 
+ 2 211 840 + 2 211 840 


+ hi —cos 127 =. 41 cos 87 
309 657600 16 588800 


21059 cos 47 


+ 1363 | 


— cos 6 2 
Ce, m) = co 4m + 4 


+4 { 960 + 192 
6 —eos10n _ 13 cos 67 11 cos 27 
+h 80640 96000 + 17280 


g f cos 127 23cos8y 1 
+8 10321920 + 6048000 92160 


10f -cos14n 53-cos 107 
+h 11857945600 + 1032192000 


4037 cos 67 439-cos 27 
+ 2419200000 + 62208000 Bu 


und die Funktionen höherer Ordnung folgen aus: 


1) E. Mathieu, a. a. O. S. 132. 


79626 240 22184) to 


bea 
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Ce,(h?,n) = cosg + + 
cos(g + 4)7 cos(g — Mn 
+ Rethg+y 32g — — 2) 
+ he { —cos (g + 6) 7 —(9? +49 + T)cos(g+2)n 
27-39 + 1)\G+2)(g+3) 
+ (9?—49+17). cos(g—2)n + cos (g — 6) 7 


cos (g + 8) 7 
+ (9° + 79? + 20g + 20) cos(g + 4) 7° 
28-3-(g + 1)°-(g— + 2): (9 + 8) 
(9? — Tg? + 20g — 20) cos(g — 4)7 
+ 
cos (g — 8)7 

Aus den angeschriebenen Reihen sieht man folgendes: 
Der Koeffizient «, in (11) ist bei Ce, gleich eins, bei Ce, 
proportional zu h?, bei Ce, proportional zu A*, bei Ce, pro- 
portional zu Ah’, usw. Er verschwindet also mit h von g-ter 
Ordnung. 

Bei den Koeffizienten «, ist immer einer, der nicht ver- 
schwindet, wenn A nach Null strebt und zwar ist es immer «,, 
wenn 2n die Ordnung der Mathieuschen Funktion bezeichnet. 
| Wir betrachten nun die Koeffizienten in der Reihe (10) 


für u. 
Wenn Ah nach Null strebt, streben sämtliche Ausdrücke 


an; (n + 0) 
nach Null. Nur der Ausdruck 


strebt nach Eins. 
Dagegen streben sämtliche Ausdrücke 
- an; 
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für 

h—>0 
nach dem Wert }. Der Ausdruck 


strebt für 
h—>0 
nach Eins. 
Diese Überlegungen werden wir im nächsten Abschnitt 
verwenden. 


6. Vergleich mit dem Kreisquerschnitt 


Heine!) hat gezeigt, daß sich die Mathieuschen Funk- 
tionen nach Besselschen Funktionen entwickeln lassen. 

Wir gehen dabei aus von der Tatsache, daß sich jede 
gerade Funktion ~(&,7) nach Mathieuschen Funktionen ent- 
wickeln läßt gemäß: 


(12) En = Ce (4,8) Coy (2). 


Als Funktion ~ nehmen wir den Ausdruck 


cos (i @ x); i=)-1, 
der ein Integral der Gleichung (1) darstellt und multiplizieren 
nun in 12 beiderseits mit Ce, (h?, n), wonach wir von — bis 
r mitteln. 
Es entsteht 


Beachtet man nun die Beziehung (11) sowie die bekannte 
Gleichung ?): 


J,,,(2) = [eos (2 cosn)- cos(2n n).dn, 


-n 


1) E. Heine, a.a. O. 8. 414. 
2) Nielsen, Handbuch der Zylinderfunkt., S. 60. 
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so geht hervor 


A, , Ge, ,,(h?, 8) 
1) - Jy (iwc + 
1 


Die Koeffizienten «, und «, sind sämtlich bekannt aus den 
Reihen des vorigen Abschnittes, Die Koeffizienten 4,,, fallen 
bei allen Anwendungen von (13) heraus. 


Wenn die Exzentrizität ce der Fundamentalellipse nach 
Null strebt, während die Konstante & einen endlichen Wert 
annimmt, so strebt 7 nach dem Polarwinkel » und die Ausdrücke 


c-Coj&E und c-Ging 


streben beide nach r, so daß die Transformationsgleichungen (2) 


übergehen in 
z=r.008p, 


y=r-sing. 
Insbesondere ist also zu beachten, daß das Produkt 
c-Coj& bzw. c- Sing 


stets endlich bleibt bei diesem Grenzübergang. Wenden wir 
diese Überlegungen, sowie jene des vorigen Abschnittes auf die 
Gleichungen (13) und (10) an, so ergibt sich: 


und 
[Ay Ge, (h?, = (iar). 


Durch Vertauschen von £ mit in erhält man: 


A,,,* Ce, (R?, m) = + (= 6-6-0081) 
und 
Ce, ,(h?, = cos 2ng, 


Ce, (h?, )e—=1. 
Bezeichnen wir mit R den Wert 


(c+ Gof ’ 
c->0 
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also den äußeren Leiterhalbmesser bei verschwindender Ex. 
zentrizität, so schreibt sich die Gleichung (10) für 


h?—->0: 
(10’) In(iar) 


um 


Es ist dies in der Tat der Ausdruck, der sich bei der 


Lösung der Gleichung (1) im Kreis mit der Randbedingung, 
daB u auf dem Kreisumfang den Wert K annehme, ergibt. 


7. Entwicklung von u bei verschwindend kleiner Exzentrizität 


Wir entwickeln vu unter Beibehaltung von A? und Ver- 
nachlässigung der höheren Potenzen. 

Zunächst ist aus früheren Überlegungen klar, daß wir in 
diesem Falle nur die ersten zwei Glieder der Reihe (10) brauchen, 
da die nächsten Glieder h? nur in höherer als der ersten 
Potenz enthalten: 


n)dn=1, 


-n 


a h* 
n)-dn= +... usw, 


während im Nenner jedes Gliedes von (10) der Ausdruck 


-an 


nicht mit A? nach Null strebt, sondern stets endlich bleibt und 
zwar gilt mit der erwähnten Näherung: 


=1, 


[Ce, (h?, + dy = usw. 


ausc 


(10” 


Fi 


Nu 


Im 
= 
|. 
= 
der 
= 
“ Sc 
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Im Nenner vernachlässigen wir additive Glieder mit A? als 


Faktor, also: 
A, + Ge, (h®, &)= R), 
A, + Ge, §,) =— J, R). 
Nach Anwendung ähnlicher Überlegungen auf die Zähler- 
ausdrücke in (10) erhält man 


J,liar).h? 
(10”) Jy (t-a+R) Bh Gea R) 
Jy (t-a+r)+h*+ cos 2p 
+ 2 b+... 


Hiermit ist die Entwicklung von x ermittelt. 


8. Das Flächenintegral von u 


Für spätere Berechnungen brauchen wir das über die 
Ellipsenfläche erstreckte Integral: 


fay-«=M. 


Für das Linienelement erhält man in elliptischen Koordinaten 
den Ausduck 


(ds)? = (Coj? — cos? + dn}. 
Das Flächenelement lautet somit: 
dz.dy = c?(Coj?& — cos? n)-d&-dn, 
und das Flächenintegral: 


M= — cos?7)- u 
0 _n 
Nun genügt aber der (4), d.h. 
Hierdurch schreibt sich das Integral: 
M= f dé = 


Schreiben wir für die Koeffizienten unter dem Summenzeichen 
in Gleichung (10) zur Abkürzung B,,, und bezeichnen wir mit 
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den oberen Strichen die Differentiation nach &, bzw. nach 1, 
so ergibt sich: 


& a 
M= = > B,, 
0 0 
(h?, &)- Ce, (h?, 1) + Ce, (h?, &) + Ce) (h?, 
Nach einem früheren Satz ist 
(02, n) = 0 


und ebenso gilt: 
Ce,,(,0)=0. 


Das Integral wird demnach: 
1) (04, &)-[ n)-dn, 
0 


oder ausgeschrieben: 


a 2 
Ce’, (h, &,) Ce, (h?,n)* an} 
0 Ce, (h?, §,) 


K 
(14 ) M = = ® z 


Wir entwickeln nun diesen Ausdruck für kleine A? unter 
Beibehaltung von h? und Vernachlässigung der höheren Po- 
tenzen. 


9. Reihenentwicklung für M 


Aus früheren Überlegungen ist sofort klar, daß wir nur 
das erste Glied der Reihe (14°) brauchen, da das zweite Glied 
bereits A* als Faktor besitzt. 

Mit dem oberen Strich ist in 14’ die Differentiation nach 
& bezeichnet worden. Nun gilt aber 


d d 
Ten Sins, 


und folglich im Grenzfall 4? —> 0: 


(1 4” 


| .. 
Der 
mit 
von 
M: 
übe 
| das 
wel 
wol 
übe 
un 
Or: 
de = a Da 
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Beachtet man dies, sowie die Überlegungen des Ab- 
schnittes 5 und 6, so erhält man 


a 


R) 
(14") (ier) h? 
{ d(ior) |r=R 5 | d(tar) ian +. 
Der Ausdruck (14”) läßt sich auch aus der Formel (10”) ableiten 
mit Hilfe von 
R 
M= |dy- |r-u-dr. 
Hierbei tritt das Integral: 
auf, für das man nach einer leichten Rechnung erhält 
Sr -- 


Unter Benutzung dieser Formel ergibt die Integration 
von (10”): 


R dJy (iar) h? [ dJ, (tar) 


K-2 


Nun ist aber die Exzentrizitätc nach Voraussetzung gegen- 
über dem Leiterradius R zu vernachlässigen und folglich ist 
das zweite Glied gegenüber dem ersten zu vernachlässigen, 
wenn wir nur Glieder mit der ersten Potenz von h? beachten 
wollen.) 

Die Integration von (10”) ergibt demnach ein mit (14”) 
übereinstimmendes Resultat. 

Hiermit sind die Berechnungen allgemeiner Art beendet 
und wenden wir uns den besonderen Fällen der Strom- 
verdrängung und der Feldverdrängung zu. 


2,2 
1) Genauer ist vorausgesetzt, daB h? = — klein von der ersten 


Ordnung gegen eins ist und e klein von der ersten Ordnung gegen R. 
Dann ist aber a® R? groß von der ersten O:duung gegen eins, g+e+d. 
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10. Stromverdrängung im elliptischen Zylinder 


Bei der Stromverdrängung verlaufen die elektrischen 
Kraftlinien parallel zur Leiterachse. Wir müssen also wu der 
elektrischen Leitungsstromdichte s gleich setzen und es ent- 
spricht der Ausdruck (14) bzw. (14”) dem gesamten elektrischen 
Strom J durch den Leiter. 

(15) M=1. 


Uns interessiert das Verhältnis der Impedanz Z pro 
Längeneinheit zum Gleichstromwiderstand R, der Längen- 
einheit. 

Bezeichnet man mit o den spezifischen Widerstand des 
Leitermateriales, so beträgt der Spannungsabfall pro Längen- 
einheit des Leiters bei Gleichstrom von der Stärke I: 


wobei a und 5 die große und die kleine Halbachse der Ellipse 
bezeichnen. 


Nach der Formel (3) gilt: 


a = c- Gof 
b=c-Ging,. 
Der Spannungsabfall bei Wechselstrom beträgt: 
e=I:-Z= 0° K. 
Das Verhältnis der Impedanz zum Gleichstromwiderstand 
wird demnach: 
Z Gof &, - Sin &, 
(16) 
wobei M der Gleichung (14’) zu entnehmen ist. 
Mit der Formel (16) ist das Problem der Stromverdrängung 
im elliptischen Zylinder formal völlig gelöst. 
Wir entwickeln nun noch Z/R, unter den gleichen Voraus- 
setzungen, die bereits früher gemacht wurden. 


Da der Nenner von M im Ausdruck Z/R, als Zähler auf- 
tritt, entwickeln wir den Nenner von M nach Potenzen von 
h? und vernachlässigen im Zähler von M bereits Glieder mit 
h? als Faktor, 


N 

V 
mit je 
i? gle 
I 
Leiter 

11. 
I 
magne 
setzen 
dissip 
2 
schla 
griff 
] 
Lang 
] 

(17) 
] 
Zylin 

(18) 
wobe; 
sieru 

(19) 
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Nach einiger Rechnung entsteht: 


2 
d(tiar) r=R 


Wie man sich leicht überzeugt, stimmt dieser Ausdruck 
mit jenem im Falle des Kreises!) völlig überein, sobald man 
h gleich Null setzt. 

Die Formel 16’ erlaubt, die Impedanz eines elliptischen 
| Leiters bei kleiner Exzentrizität zu berechnen. 


11. Magnetische Feldverdrängung im elliptischen Zylinder 


Im Falle der magnetischen Feldverdrängung verlaufen die 
magnetischen Kraftlinien parallel zur Leiterachse und wir 
setzen die Größe u der magnetischen Feldstärke H gleich. 

Uns interessiert in diesem Falle die in dem Zylinder 
dissipierte Joulesche Energie, sowie die in ihm enthaltene 
magnetische Energie. 

Zur Berechnung dieser Größen kann man drei Wege ein- 
schlagen.?) Wir benutzen jene Methode, welche auf dem Be- 
griff der induzierten elektromotorischen Kraft beruht. 

Die Magnetisierungsspule besitze w Windungen in der 
Längeneinheit und führe den effektiven Strom I. 

Dann gilt: 

(11) gu 


Der gesamte magnetische Fluß durch den elliptischen 
Zylinder beträgt: 
18) F=u-M, 
wobei M der Gleichung (14°) zu entnehmen ist. 

Dieser Fluß induziert in der Längeneinheit der Magneti- 
sierungsspule eine elektromotorische Kraft: 

1 OF P w 

Man muß nun bedenken, daß F letzten Endes in J aus-. 
gedrückt ist, vermöge der Gleichungen (18) und (17). 


1) E. Jahnke u. F. Emde, Funktiontafeln, 8. 143, Gleichung 9. 
2) M. Strutt, Ann. d. Phys. 82. 8. 605. 1927. 
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Der reelle Teil des Ausdrucks (19) stellt eine Spannung 
dar, die genau entgegengesetzte Phase besitzt, wie der Strom I. 
Dieser reelle Teil, multipliziert mit J, gibt also die in dem 
elliptischen Zylinder pro Länßeneinheit dissipierte Joulesche 
Energie an. 

Man findet somit für die Joulesche Wärme: 

a 2 


(20) Q= 0-w?- I?.«?- Re £ Ge, 


’ 


f [Ce, n))’dn 


wobei Re( )=reeller Teil. 


Verwendet man für o und J die praktische Einheit und 
für die Längeneinheit das cm, so erhält man aus der Glei- 
chung (20) sofort Q in Watt/cm. 

Bei kleiner Exzentrizität können wir die Reihe nach # 
entwickeln und es entsteht gemäß der Gleichung (14”): 


= 0:w2.I2. (dd Ger) 
ao) Q = 0-w?- I Re {i 


Der imaginäre Teil des Ausdrucks (19) stellt eine Spannung 
dar, deren Phase um 90° gegen diejenige des Stromes J ver- 
schoben ist. Dieser Teil, multipliziert mit Z, gibt demnach die 
in der Längeneinheit des elliptischen Zylinders enthaltene 
magnetische Energie an: 


= 2 
(21) T= Ce’, » (h®, &) | (hi nar} | 


GE) 
0 S [Cen (h*, dy 


wobei Jm( ) = imaginärer Teil. 


Ähnlich wie für die Joulesche Energie erhalten wir für 
kleine Exzentrizität den Ausdruck: 


J, (ie R) diier) 


_ # 

Hiermit ist das Problem der magnetischen Feldverdrängung 
im elliptischen Zylinder gelöst. 
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12. Schlußbemerkungen 


Eigenschaften der Mathieuschen Funktionen, deren wir 
uns zur Lösung der vorliegenden Randwertaufgabe bedienten, 
sind nur wenige bekannt. Außer den Orthogonalitätsrelationen 
lassen sich noch einige Eigenschaften unmittelbar aus den 
Reihenentwicklungen ablesen. Wir gaben sie in den Ab- 
schnitten III, V, VI und VII an und benutzten sie zugleich. 

Mit Wirbelströmen im elliptischen Zylinder beschäftigte sich 
bereits R. Maclaurin.!) Seine Betrachtungen waren jedoch 
auf die Auffindung der Eigenzeitkonstanten solcher Zylinder ge- 
richtet, ein Problem, dessen Behandlung und Lösung vollständig 
parallel verläuft mit der Bestimmung der Eigenschwingungen 
einer elliptischen Membran?) und eines elliptischen Wasser- 
beckens.°) 

Wir möchten darauf hinweisen, daß der Ausdruck (10) auch 
die Durchbiegung einer hydrostatisch belasteten Membran angibt, 
die, wie wir an anderer Stelle*) zeigten, in einfacher Be- 
ziehung steht zu den in dieser Arbeit behandelten Problemen. 

Für die numerische Beherrschung der Mathieuschen 
Funktionen ist zurzeit noch wenig geschehen. 

Zunächst sind ihre Reihenentwicklungen nur für kleine A 
bekannt. Von keiner der Reihen ist das allgemeine Glied be- 
kannt (obwohl man dies bei obertlächlicher Betrachtung von 
Whittakers°®) Reihen vermuten könnte), Eine Tabellierung 
ist bis jetzt sogar für diese kleine A nicht erfolgt. Für die 
hier behandelten Probleme wäre eine Tabellierung der 
Mathieuschen Funktionen im engeren Sinn für imaginäre Ah? 
notwendig. 

Unsere Entwicklungen (10”, 14”, 16, 20’) und (21’) sind 
numerisch nicht ohne weiteres verwertbar, da es keine Tabellen 
der Besselschen Funktion J,’ mit komplexem Argument gibt. 
Man kann aber mit Hilfe der Formel: 
dJ,(i-a@-R) 4 dJ, (ia R) _2i R) — ad J, (ia R) 
d(i-a-+ R) =~ RP d(ta R) aR o (ta R) d(ia R) 


1) R. Maclaurin, Proc. Camb. Phil. Soc. 17. S. 41. 1899. 
2) E. Mathieu, a.a. O. 
3) H. Jeffreys, Proc. Lond. math. Soc. 23. S. 440. 1925. 
4) M. Strutt, Physica 6. S. 364. 1926. 
5) E. T. Whittaker, a.a.0. 8.411. 

Annalen der Physik. IV. Folge. 84. 33 
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diese Funktionen durch die Funktionen J, und J,’ ausdrücken, 
die für komplexes Argument x-(1 + i) tabelliert sind.') 

Für sehr große Werte 

ak 

lassen sich an Stelle von (14”, 16’, 20’, 21’) Näherungsformeln 
schreiben, die man erhält, wenn man die Besselschen Funk- 3 
tionen durch ihre asymptotischen Werte fiir groBes Argument 
ersetzt. 


Eindhoven, Natuurkundig Laboratorium der N. V. Philips 
Gloeilampenfabrieken, August 1927. 


1) H.B. Dwight, Journ. Am. Inst. El. Eng. 42. $. 830. 1923. ges 


(Eingegangen 12. September 1927) 
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3. Über Koppelschwingungen kontinuierlicher 
Teilsysteme'); 
von E. Waetzmann und K. Schuster 


$ 1. Fragestellung 

Die Untersuchung von Schwingungen mechanischer Systeme 
gestaltet sich — namentlich in theoretischer Hinsicht — ver- 
schieden, je nachdem, ob es sich um diskrete Massenpunkte 
oder um Kontinua handelt. Mathematisch kommt der Unter- 
schied dadurch zum Ausdruck, daß man es im ersten Falle 
mit gewöhnlichen, im zweiten mit partiellen Differential- 
gleichungen zu tun hat. Die Bestimmung der Eigentöne, die 
das Hauptziel jeder Schwingungsuntersuchung bildet, führt dort 
auf eine algebraische, hier auf eine transzendente Frequenz- 
gleichung. Außer den Eigentönen spielt bei den Schwingungen 
der Kontinua noch die räumliche Amplitudenverteilung, die 
einem bestimmten Eigentone zugeordnet ist, eine Rolle: die 
Eigenfunktion. Sie äußert sich in Erscheinungen wie den 
Chladnischen und Kundtschen Figuren. 

Ist ein Schwingungssystem aus zwei oder mehreren Teil- 
systemen zusammengesetzt, die wenigstens angenähert auf- 
einander abgestimmt sind, so spricht man von „gekoppelten 
Schwingungen“ oder von „Koppelschwingungen“. Ihre Theorie 
ist für den speziellen Fall zweier Massenpunkte mit je einem 
Freiheitsgrade von M. Wien?) entwickelt worden. Ent- 
sprechende Messungen, die gute Übereinstimmung mit der 
Theorie ergaben, wurden am mehrfachen Pendel ausgeführt.) 
So hatte man seit der grundlegenden Arbeit M. Wiens über 


1) Vorgetragen in der Sitzung des Gauvereins Thüringen-Sachsen- 
Schlesien der D. Phys. Ges. am 12. Juni 1927. 

2) M. Wien, Wied. Ann. 61. $. 151, 1897. 

3) E. Barton u. H. Browning, Phil. Mag. (6) 34. S. 246. 1917; 
40. 8.611. 1920; 42. 8.17. 1921; A. Narayan, a. a. O. 43. S. 567. 
1922; C. Lees, a. a. O. 48. S. 129. 1924; V. Solovieff, a.a.O. 50. 
8. 612. 1925. 
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die Koppelschwingungen von Gebilden, die als Massenpunkte 
behandelt werden dürfen, einigermaßen klare Vorstellungen, 
Allerdings kommt man leicht in Verlegenheit, wenn man in 
speziellen Fällen etwa angeben soll, ob Kraft- oder Massen- 
koppelung vorliegt, und wie die Koppelungsverhältnisse und 
ihre Wirkungen im einzelnen beschaffen sind. 

Es liegen auch schon einige experimentelle Untersuch- 
ungen über die Schwingungen gekoppelter Kontinua vor.') 
Zur Berechnung benutzte man dabei die Ergebnisse der Massen- 
punktstheorie. Das genügt in manchen Fällen, aber nicht 
immer. Die typischen Koppelungserscheinungen, vor allem die 
veränderte Lage der Eigentöne des Gesamtsystems gegenüber 
denjenigen der Einzelsysteme, machen sich auch bei Schwin- 
gungen gekoppelter Kontinua bemerkbar. Aber man kann nicht 
erwarten, daß die an Hand der Massenpunktsschwingungen 
hergeleiteten Begriffe und Gesetze im einzelnen auf den kom- 
plizierteren Fall der Kontinuumsschwingungen übertragen werden 
können. Ein wirklich klarer Einblick in die Vorgänge ist 
durch die rohe punktförmige Betrachtungsweise nicht zu ge- 
winnen. Deshalb haben wir ein paar auch praktisch wichtige 
Modelle gekoppelter Teilsysteme, an denen sich der Übergang 
vom punktförmigen zum kontinuierlichen Verhalten gut ver- 
folgen läßt, systematisch untersucht. In der vorliegenden 
Arbeit berichten wir über die experimentellen Ergebnisse und 
ihre Deutung, während die theoretische Berechnung in einer 
gesonderten Arbeit von K. Schuster gegeben wird. 

Als Teilsysteme wurden benutzt: erstens Telephon- 
membranen (TM.), die den Vorteil bieten, daß eine elektrische 
Meßmethode zur Bestimmung der Eigenschwingungszahlen an- 
gewandt werden kann, und zweitens zylinderförmige Lufträume 
mit starren Seitenwänden. Die Höhe A der Luftsäule wurde 
von 1 cm bis zu einer halben Luftwellenlänge 2 des Grundtones 
der TM. variiert. Folgende Systeme wurden untersucht: 
1. eine Luftsäule von überall gleichem Durchmesser (5,3 cm), 
die an dem einen Ende durch eine TM., am andern durch 
eine starre Wand verschlossen ist (System I, Fig. 1a), 2. eine 


1) H. Vogel, Ann. d. Phys. 62. S. 247. 1920; W. Hahnemann, 
H. Hecht und H. Lichte, Ztschr. f. techn. Phys. 4. S. 93. 1923. 
E. Mallet und G. Dutton, Journ. Inst. Electr. Eng. 63. S. 502. 1925. 
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Luftsäule von kleinerem Durchmesser (1 cm), die einerseits 
wieder durch eine starre Wand verschlossen ist, und am 
andern Ende in eine Kapsel größeren Querschnittes mündet, 
die ihrerseits durch die TM. abgeschlossen 
ist. Die Höhe der Kapsel (1,5 cm) wurde 
bisher nicht variiert (System II, Fig. 1b), a 
3. eine Luftsäule von überall gleichem 
Durchmesser (5,3 cm), die an beiden +” 


Enden durch TM. verschlossen ist, die 
ungefähr den gleichen Grundton haben [) 
(System III, Fig. 1c). h 


$ 2. Meßmethode und Apparatur 


Eines der Telephone wurde durch e 
einen stetig veränderlichen Röhren- 
generator (Rohrsummer nach Siemens & - A ” 
Halske) erregt, und es wurde die Fig. 1 
Summe der im Eisen verbrauchten 
Leistung Z, und der mechanisch-akustischen Leistung Z bei 
verschiedenen Frequenzen nach Hahnemann und Hecht!) be- 
stimmt, womit die Eigenschwingungszahlen des Systems gegeben 
sind. Sie fallen mit den bei konstanter Amplitude der er- 
regenden Kraft bestimmten Resonanzstellen der mechanisch- 
akustischen Leistung fast genau zusammen. 

Die Summe (Z, + Z) steht mit der Kreisfrequenz w in der 
Beziehung o(W - Wo) 

(! 


worin Wg den Gleichstromwiderstand, W = W(w) den Watt- 
widerstand und / = /(w) die Selbstinduktion des Telephons be- 
deuten. ce wird für die Berechnung als konstant angesetzt, 
was darauf hinauskommt, daß die Amplitude der erregenden 
Kraft konstant gehalten wird. Da es sich nur um Relativ- 
messungen handelt, setzen wir c=1. Zur Messung der 
Frequenz diente ein Frequenzmesser nach Siemens & Halske. 
Neben dem Gleichstromwiderstand müssen dann noch W und / 
gemessen werden, um die (Z,+ Z)— »-Kurven zu erhalten. 
Da die Eisenverluste in erster Näherung der Frequenz pro- 


1) W. Hahnemann u. H. Hecht, Ann. d. Phys. 60. S. 454. 1919. 
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portional sind, erhält man Kurven von der Gestalt der in Fig. 2 5,3 
dargestellten Kurve. Die der geraden Linie zugehörenden den 
Ordinatenwerte sind von den Gesamtwerten der Ordinaten ab- 
zuziehen, um die Resonanz- 
kurven für die mechanisch- 
akustische Leistung allein 
zu bekommen. 

Fig. 3 zeigt das Schal- 
tungsschema für die Mes- 
sung von W und. R be- 
deutet den Rohrsummer, 
F den Frequenzmesser, 
S das Schwingungssystem 
mit dem zu messenden 
Telephon (W, 1), C einen Glimmerkondensator, w,, w, und w, 
sind induktionsfreie Widerstände, 7, und 7, die Anschluß- 
stellen für das als Nullinstrument dienende Telephon. 7, ist 
die Anschlußstelle in der Meßbrücke, 7, die für den Wagner- 
schen Hilfszweig!) HH mit der Erdung Z. Dieser Hilfszweig 
wurde zur Beseitigung kapazitiver Störungen angewandt. Um 
induktive Störungen zu vermeiden, wurden die Endpunkte der 
Brückenanordnung dicht aneinandergelegt?); die Widerstände 
selbst standen möglichst weit voneinander entfernt; die Zu- 
leitungsdrähte waren bifilar gewickelt. Der Brückenstrom ver- 
schwindet, wenn 


w 
Fig. 2 22 


(2) W-w, 


ist, wobei alle Größen in CGS,-Einheiten auszudrücken sind. 
Daraus ergeben sich die gesuchten Größen # und !: 
(3) W= l= w,-w,-C. 

Die Verbindung der einzelnen Teile der zu untersuchenden 
Systeme mußte so gestaltet werden, daß die Telephone selbst, 
inbesondere die Einspannung der TM., beim Anschließen an 
den Luftzylinder möglichst keine Veränderung erfuhren. Hierzu 
wurde folgende Anordnung (Fig. 4) getroffen: Die TM. von 


1) K. W. Wagner, Elektrot. Ztschr. 32. S. 1001. 1911. 
2) E. Griineisen u. E. Giebe, Wiss. Abh. d. Phys.-Techn. Reichs- 
anstalt 5. S. 83, 1921. 
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5,3 cm freiem Durchmesser und 0,01 cm Dicke ist zwischen 
den beiden Messingringen A und B eingeklemmt. Der innere 
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Ring 4 ist auf das Telephongehäuse aufgeschraubt, der äußere 
B ist flanschartig erweitert und kann auf den Messingring C, 
aufgeschraubt werden, in den das Rohr Z endigt. Dieses be- 
stand bei kleiner Höhe k (bis zu 5 cm) aus einem festen Stück, 
bei größeren Höhen wurden zwei ineinander verschiebbare 
Rohre verwendet. Auf den Endring C, des Messingzylinders 
wurde bei System I eine starre Platte P aufgeschraubt; beim 
System II eine starrre Platte, in die das dünne auf der andern 
Seite verschlossene Messingrohr Q endet; beim System III ein 
zweites Telephon. ; 

Die Frequenz konnte auf etwa 2 Hertz genau gemessen 
werden. Ebenso genau konnten die Maxima der mechanisch- 
akustischen Leistungen, also die Eigenschwingungszahlen, aus 
den Kurven abgelesen werden. Die Eigentöne der Telephone 
blieben aber bei den Messungen nicht genau konstant, sondern 
änderten sich (auch bei möglichst gleichmäßigem, nicht zu 
festem Anziehen der Schrauben) immer noch in unsystema- 
tischer Weise um Beträge von etwa 5 Hertz. Gelegentlich 
kamen auch bedeutend größere Sprünge vor, die sich natür- 
lich sofort bei der Messung bemerkbar machten, und gelegent- 
lich zeigte sich viele Stunden hindurch gute Konstanz. Damit 
ist die Genauigkeit der Meßresultate umschrieben. 


§ 3. Meßergebnisse 


In den Figg. 5—10 sind einige Proben der (L,+ L)—o- 
Kurven zusammengestellt. Als Abszissen sind nicht die Kreis- 
frequenzen, sondern die Frequenzen in Hertz aufgetragen. Der 
Frequenzbereich ist so gewählt, daß sich im allgemeinen nur 
die Grundtöne der Teilsysteme bzw. die ihnen entsprechenden 
Eigenténe des Gesamtsystems bemerkbar machen können. 
Andernfalls würden sich Schwierigkeiten ergeben, weil für die 
Berechnung der Eigentöne Frequenzunabhängigkeit des Mem- 
branprofils angenommen ist. 

Für die Deutung der Resultate ist es wichtig, auch die 
Eigentöne, namentlich die Grundtöne der „bloßen Telephone“ 
zu kennen, d.h. die Grundtöne der mit ihren Gehäusen ver- 
bundenen TM., wenn von der äußeren Seite der Membranen 
her gar keine Kräfte wirken. Dieser Zustand läßt sich streng 
nicht verwirklichen. Denn wenn man das Telephon im Vakuum 


beoba 
norm: 
freien 
Luftn 
Hahı 
Luftn 
(4) 
wobei 
der ( 
zität 
[7 
die 
Eige 
6) 
Hier 
Mas: 
die 
R= 
der 
tieft 
töne 
lich 


Über Koppelschwingungen kontinuierlicher Teilsysteme 513 


beobachten wollte, so würde der im Gehäuse herrschende 
normale Luftdruck die Membran nach außen durchbiegen. Im 
freien, unbegrenzten Luftraume bewirkt die mitschwingende 
Luftmasse eine Vertiefung des Grundtones, die sich aber nach 
Hahnemann und Hecht!) berechnen läßt. Die mitschwingende 
Luftmasse ist 

(4) Am=0,4- 0, 

wobei R der Radius der TM. und o die Luftdichte ist. Da 


der Grundton des Telephons w = V — ist, worin c die Elasti- 
ztät und m die Schwingungsmasse der Membran bedeuten, ist 


1 N 4 
600 800 900. 1000 1200 7300 *0OHerrz 
Fig. 5 


die dem Massenzuwachs 4m entsprechende Anderung des 
Eigentones 
dm 1 0,4-R% Re 


1 


Hierbei ist das Verhältnis von Schwingungsmasse zu Gesamt- 
masse mit Hahnemann und Hecht gleich } gesetzt. « ist 
die Membrandichte, D ihre Dicke. In unserem Falle, wo 
R= 2,65 cm und D = 0,086 gem”? war und der Grundton 
der Telephone bei ungefähr 700 Hertz lag, beträgt die Ver- 
tiefung also etwa 9 Hertz. Um diesen Betrag sind die Grund- 
tine, die aus den Leistungskurven der im Luftraum befind- 
lichen Telephone gewonnen waren, zu erhöhen. 


1) W. Hahnemann u. H. Hecht, Phys. Ztschr. 18. 8. 261. 1917. 
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Zur angenäherten Bestimmung der Obertöne wurde Kork- 
staub auf die TM. gestreut, und es wurde das Auftreten der 
Staubfiguren beobachtet. Der Oberton mit einer Durchmesser- 
knotenlinie wurde überhaupt nicht festgestellt, was bei der 


800 900 Hertz 
Fig. 8 
| 
700 800 Hertz 
Fig. 9 


Konstellation Telephonmagnet-Membran wohl nicht wunderbar 
ist. Derjenige mit zwei Durchmesserknotenlinien war ziemlich 
schwach und lag in der Gegend von 1250 Hertz. In den 
Leistungskurven wurde seine Einwirkung nicht beobachtet, was 
damit zusammenhängen dürfte, daß sich hier ein Maximum 
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nur dann bemerkbar macht, wenn die mechanisch-akustische 
Leistung eine gewisse Größe im Verhältnis zu der im Eisen 
verbrauchten Leistung erreicht. Sehr stark war dann der 
dritte Oberton mit einer Kreisknotenlinie, der erste zentral. | dem 
symmetrische Oberton, der in der Gegend von 1600 Hertz lag, | dopp 
Er ergab auch in der Leistungskurve des Telephons bei Er- | zwecl 
regung durch 1600 Hertz ein starkes Maximum. Aber auch § für e 
bei Erregung in der Grundfrequenz 800 Hertz scheint er zu 
gelegentlichen Störungen Veranlassung gegeben zu haben. Der § noch 
erregende Wechselstrom war so weit sinusförmig, daß in den schw: 
Leistungskurven der Telephone für sich bei Erregung in der stelle 
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genannten Grundfrequenz keine Andeutung eines Maximums 
gefunden wurde. Sobald aber bei einem unserer Gesamt- 


d 
systeme I, II oder III ein Eigenton in der Nähe dieser Stelle at 
lag, wies die Resonanzkurve eine Störung auf. Besonders “_~ 


deutlich zeigt sich dieser Effekt in der dritten Kurve der ’ 
Fig. 6. Hier lag der Grundton der TM. bei etwa 705 Hertz; i 
ferner ist zu beachten, daß bei der Koppelung auch der Ober- Rig 
ton eine leichte Verschiebung erfahrt. TM. 

Die Figg. 5 bzw. 6 stellen die MeBresultate nur für einige Wi 
größere h-Werte an den Systemen I bzw. II dar. Die Figg. 1, für 
8 und 9 zeigen die Resonanzkurven des Systems III für die 
Höhen 1, 2 und 3cm. Von den beiden Kurven einer Figur 
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bezieht sich die eine auf Primärerregung des einen, die andere 
auf Primärerregung des zweiten Telephons. Die Zuleitung zu 
dem nicht primär erregten Telephon war immer geöffnet. Die 
doppelten Kurvenaufnahmen sind hauptsächlich zu Kontroll- 
zwecken gemacht worden. Fig. 10 enthält die Meßresultate 
fir einige größere h-Werte am System III. 

Zur bequemeren Diskussion der Ergebnisse wollen wir 
noch für jedes System den Zusammenhang zwischen den Eigen- 
schwingungszahlen = und den A-Werten in Kurvenform dar- 
stellen. Bevor wir die Diskussion beginnen, miissen wir aber 
noch kurz auf die theoretische Behandlung des Problems 
eingehen. 

§ 4. Theoretisches ') 

Wir denken zunächst speziell an das System III. Es 
wird für jede TM. eine Kraftgleichung aufgestellt mit den 
elastischen Kräften und den von den Schwingungen der Luft- 
siule herrührenden Druckkräften. Das Profil der TM. wird 
von vornherein in bestimmter Form angesetzt, so daß die 
Druckkräfte durch die Mittelpunktselongationen der Membranen 
ausgedrückt werden können. Auf diese Weise nimmt die 
Gleichung für jede TM. die Gestalt einer Wienschen 
Koppelungsgleichung an. Das bringt allerdings im allgemeinen 
Falle höchstens einen formalen Nutzen, weil die Koeffizienten 
der einzelnen Glieder der Gleichungen — im Gegensatz zu 
Wien — noch frequenzabhängig sind, und damit eine an- 
schauliche Deutung unmöglich machen. 

Nur für sehr kleine h-Werte(k < A) werden die Koeffizienten 
der für unser System III aufgestellten Gleichungen frequenz- 
unabhängig. Dann läßt sich die Frequenzgleichung leicht nach 
®, auflösen und die Bestimmungsstücke für w, lassen sich in 
einfacher Weise ausdrücken durch die Eigentöne %,/2r und %,/2 x 
der beiden TM. und — was noch interessanter ist — die 
Eigentöne von erster TM. + Luftsäule %,/2” und von zweiter 
TM. + Luftsäule 7,/2%. Es folgt hieraus in Analogie zu 
Wienschen Gleichungen und Betrachtungen, daß das System III 
für kleine h-Werte als gekoppelt aus den beiden Teilsystemen 


1) Über alle Einzelheiten vgl. die nachstehende Arbeit von K. 
Schuster. 
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erste TM. + Luftpolster (nicht TM. für sich) und zweite TM. + Luft. 
polster anzusehen ist. Übrigens läßt sich dieses Resultat 
auch ohne jede Theorie aus unseren Meßergebnissen ablesen. 
Die — wiederum in Analogie zu Wienschen Rechnungen — 
sich ergebende Formel für die Eigenschwingungszahlen lautet: 


N.? 9.2 

(6) ofa + 

Bevor wir auf den Fall auch größerer h-Werte und damit 
auf die eigentliche Kontinuumstheorie eingehen, wollen wir 
einige unserer Messungen für kleine h-Werte an Hand der 
Formel (6) prüfen. Die Werte %,, %,, 2, und 7, wurden aus 
dem Leistungsverbrauch gemessen, desgleichen die beiden 
Eigenschwingungszahlen », des aus den beiden Telephonen mit 
den Grundtönen %, und %, kombinierten Systems III. Außerdem 
wurden die w,-Werte für System III nach Formel (6) berechnet. 
In der Tabelle sind die Werte für drei verschiedene Höhen 
zusammengestellt. Die Messungen für jede Einzelhöhe wurden 
zu ganz verschiedenen Zeiten und unter teilweise verschiedenen 
Bedingungen ausgeführt. 


h Ne n, Ng @, ber. 

1 688 720 843 868 700, 979 703, 985 

2 688 696 781 7171 684, 853 692, 850 

3 686 700 156 764 | 696, 809 | 691, $25 


Die Übereinstimmung zwischen beobachteten und be- 
rechneten Werten ist für unsere Zwecke durchaus befriedigend. 

Für größere Ah-Werte erhält man als Bestimmungs- 
gleichungen für ®, beim System I 


beim System II 
2 
(8) = — 
r 2 (m +2) 
Mm, 


beim System III 


(9) (cotg hh + 
Hierin bedeuten w, die Eigenfrequenzen der Gesamtsysteme, 
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n den Eigenton der bloßen TM., o die Luftdichte, a die Schall- 


geschwindigkeit, & ist gleich * = a e und D sind die 


a 

Dichte bzw. Dicke der TM., a, bedeutet das Durchbiegungs- 
volumen der TM. im Verhältnis zum Volumen des Zylinders, 
der den gleichen Durchmesser und die Mittelpunktsamplitude 
der TM. zur Höhe hat, u den Massenfaktor, das ist das Ver- 
hältnis von Schwingungsmasse zu Gesamtmasse, A, den Radius 
des engen Rohres, R den der Kapsel, h, ihre Höhe, k, die von 
ihr gelieferte Zusatzelastizität, und y ist eine Größe, die von 
der Art der Luftschwingungen in der Grenzebene von Kapsel 
und engem Rohr abhängt. 

Wenn wir wie geschehen das Profil der TM. von vorn- 
herein festlegen, so begehen wir einen gewissen Fehler und 
sind auf kleine Frequenzbereiche beschränkt. Man kann aber 
auch korrekter vorgehen. Es werden zunächst nur gewisse 
Annahmen über den Zusammenhang der Deformationen der 
TM.und der elastischen Kräfte gemacht, während eine Änderung 
des Profils mit der Frequenz und durch den Einfluß der Luft- 
schwingungen noch nicht ausgeschlossen wird. Für die Luft- 
siule wird die Wellengleichung aufgestellt, und die Telephon- 
gleichungen werden durch die Bedingung, daß an den Ober- 
flächen die Elongationen von TM. und Luft die gleichen sind, 
zu Randbedingungen der Wellengleichung. Die Berechnung 
der Koppelschwingungen wird mit diesem Ansatz zu einer 
Randwertaufgabe, deren Lösung nach der Ritzschen Methode 
möglich ist und die Vorgänge genau klärt. Beschränkung auf 
die erste Näherung in der Ritzschen Lösung bedeutet, daß 
das Profil noch als frequenzunabhängig und als unabhängig 
von dem Einfluß der Luftsäule angenommen wird. In der 
Tat ergibt die Durchrechnung für die Systeme I und III 
Frequenzformeln, die mit den Formeln (7) und {9) bis auf den 
Faktor a,?/u übereinstimmen. Für diesen Faktor, dessen Be- 
trag wir nach Hahnemann und Hecht zu 5/, angesetzt haben, 
tritt hier der Faktor '/, ein. 


$5. Diskussion der Ergebnisse 


In den Kurven der Figg. 11—13 sind für die Systeme I, 
II und III die nach den Formeln (7, 8 und 9) berechneten 
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o,-Werte als Funktionen der Höhen der Luftsäulen zusammen. | jhn | 
gestellt. Als Abszissen sind nicht die Kreisfrequenzen, sondern | Lufts 
die Frequenzen (Hertz) aufgetragen. Die aus den Kurven der f zwei 
Figg. 5—9 und aus einigen nicht veröffentlichten Kurven ent- | sich ı 
nommenen experimentellen Werte sind als Kreise eingezeichnet. | mitei 
Dabei ist für die Kurven des Systems II noch zu bemerken, 
daß infolge einer Unbestimmtheit von y der untere Eigenton 
für h = 12 cm willkürlich auf die Kurve gelegt wird, wodurch | dem 
die Werte bei A = 10 cm stärker streuen, als nach der In. | same 
konstanz der Telephongrundtöne zu erwarten ist. beka 


Beim System I (Fig. 11) liegen die Eigentöne für kleine | [uf 
h-Werte sehr weit auseinander und sehr hoch. Der tiefere 


Is 
Eigenton ist der von TM. + Luftpolster, der höhere der des ‚oo 
Luftpolsters für sich. Bei dem großen Abstande beider findet | der 


eine gegenseitige Beeinflussung kaum statt; und die anschau- 
lichste Vorstellung ist die, daß das Ganze noch ein einheitliches Eig 
Gebilde ist, bestehend aus einer TM., deren Elastizität durch 
die Elastizität des Luftpolsters vergrößert ist und das noch | die 
einen hohen Oberton besitzt, weil das Luftpolster auch für sich | met 
schwingen kann. Tei 

Für größere Werte von A nähert sich aber der tiefere 
| Eigenton des Systems mehr und mehr dem der bloßen TM. - 
ohne Luftpolster, der bei 700 liegt, und sinkt sogar noch unter 
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ihn hinab; auch der höhere weicht von dem Eigenton der 
Luftsäule mehr und mehr ab. Jetzt müssen wir also von 
zwei gekoppelten Teilsystemen sprechen, der bloßen TM. für 
sich und der Luftsäule für sich, die durch Teile ihrer Elastizität 
miteinander gekoppelt sind. Betrachten wir speziell den Wert 


h= + = 24cm. Hier ist der Eigenton der Luftsäule gleich 


dem der TM., und durch die Koppelung wird dieser gemein- 
same Eigenton des Teilsystems aufgespalten in zwei, die in 
bekannter Weise symmetrisch zu ihm liegen. Besonders 


interessant ist noch der Fall h = = 12cm. Hier ist der 


Grundton der Luftsäule gleich 1400 Hertz, der der TM. ist 
700 Hertz, und mit diesen beiden Eigentönen der Teilsysteme 
stimmen die des Gesamtsystems so gut wie vollständig überein. 
Die beiden Teilsysteme sind entkoppelt. Diese Tatsache wird 
verständlich, wenn man sich klarmacht, daß für den Ton 700 
an dem der TM. zugewandten Ende der Luftsäule sich ein 
Schwingungsbauch ausbildet, so daß keine Druckkräfte auf die 
Membran ausgeübt werden. Für den Ton 1400 wird die gegen- 
seitige Beeinflussung der Schwingungen von Luftsäule und 
Membran deshalb verschwindend klein werden, weil der Eigenton 
der Membran zu weit abliegt. 

Beim System II (Fig. 12) greifen wir nur den Fall heraus, 
daß die Höhe A der engen Luftsäule = 10 cm ist. Hier ist 
das Gesamtsystem aufzufassen als gekoppelt aus folgenden 
zwei Teilsystemen: TM. + Luftpolster der Kapsel und enge 
Luftsäule mit einer Art Mündungskorrektion infolge des Luft- 
polsters. Das Luftpolster ist das beiden Teilsystemen gemein- 
same Koppelungselement und bewirkt außerdem, daß das Ende 
der engen Luftsäule, das der TM. zugewandt ist, nicht als 
geschlossen, sondern als offen aufzufassen ist. Dann liegt der 
Eigenton dieser Luftsäule bei 10 cm Höhe etwa bei 800 Hertz 
und das ist auch der Eigenton von TM. + Luftpolster. Und 
die beiden Eigentöne des Gesamtsystems liegen wieder sym- 
metrisch zu dem gemeinsamen Eigenton der eben definierten 
Teilsysteme; also ganz anders als beim System 1, wo wir für 


h= + Entkoppelung hatten. 


Beim System III (Fig. 13) haben wir für kleine Werte von 
Annalen der Physik. IV. Folge. 84. 34 
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h, wie im vorigen Paragraphen schon erwähnt wurde, zwei 
Teilsysteme anzunehmen: erste Membran + Luftpolster und 
zweite Membran + Luftpolster. Der gemeinsame Eigenton 
dieser Teilsysteme liegt z.B. für k = 1 cm bei etwa 850 Hertz; 


em 


1700 7200 7200 


Fig. 13 


er wird durch die Koppelung aufgespalten in die beiden sym- 
metrisch dazu liegenden etwa 700 und etwa 1000 Hertz. Das 
Koppelungselement selbst, nämlich das Luftpolster, hat für sich 
genommen noch keinen Einfluß, denn sein Eigenton liegt weit 
draußen bei etwa 17000 Hertz. 

Für größere Werte von h wird die Luftsäule aber wieder 
selbständig, und wir haben mit drei Teilsystemen zu operieren: 
erste Membran (nicht mehr erste Membran + Luftpolster, 
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sondern erste Membran für sich), zweite Membran für sich 
und Luftsäule für sich. Fir 1 = — = 24 cm fällt der Eigenton 
der Luftsäule mit dem der TM., 
etwa 700 Hertz, zusammen. Das Ge- 

samtsystem hat als Eigentöne diesen Pi 
gemeinsamen Kigenton der drei Teil- ” 
systeme und zwei symmetrisch dazu ar 

liegende. Für h = + = 12cm ist 
die Koppelung äußerst lose. Die Be- 
rechnung der Amplitudenverteilung 
in der Luftsäule ergibt für den höheren 
der zwei letztgenannten Töne einen b 
Zustand gemäß Fig. 14a und für den 

tieferen einen Zustand gemäß Fig.14b. 

In beiden Fällen grenzen also an Fig. 14 

die TM. Schwingungszustände in der 


Luftsäule, die nur um 2/8 von dem Schwingungsbauch ent- 
fernt sind, so daß keine starke Beeinflussung zu erwarten ist. 


§ 6. Zusammenfassung 


Um ein klares Bild von den Koppelungsvorgängen zu be- 
kommen, ist in erster Linie notwendig, das System bei der 
Betrachtung passend zu unterteilen. Dabei muß als Grund- 
satz gelten, daß von den einzelnen Elementen des Gesamt- 
systems nur diejenigen als selbständige Teilsysteme eingeführt 
werden, deren Eigenperioden nicht zu weit auseinanderliegen. 
Elemente mit weitabliegender Eigenperiode sind also nur als 
Teile von Teilsystemen zu betrachten. Für kleine h-Werte 
wird die Luftsäule einfach zu den TM. hinzugeschlagen, indem 
deren Elastizität um die Elastizität der Luftsäule vergrößert 
wird. Das Teilsystem Membran + Luftpolster kann dann auch 
noch als punktförmig betrachtet werden. Für größere Höhen 
ist das aber nicht mehr möglich, und der tiefere Grund ist 
natürlich der, daß die einzelnen Teilchen der Luftsäule jetzt 
ganz verschiedene Schwingungen ausführen und daß die Schwin- 
gungen ein und desselben Teilchens ganz andere werden, wenn 
sich die Frequenz irgendwie wesentlich ändert. Dann muß 
die Luftsäule als selbständiges und zwar nicht mehr punkt- 
34* 
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förmiges, sondern kontinuierliches Teilsystem eingeführt werden. 
Bei der TM. wird die Schwingung irgendeines ihrer Massen- 
teilchen erst dann eine ganz andere, wenn die Frequenz aus 
dem Bereiche eines der Telephoneigentöne bis in den Bereich 
eines anderen Eigentones schwankt. Der Frequenzbereich, 
innerhalb dessen eine punktförmige Betrachtung ausreichend 
ist, ist hier also größer als bei den Lufträumen. 

Koppelungsbetrachtungen, wie wir sie von punktförmigen 
Systemen her gewohnt ‘sind und die dort eingeführten Begriffe 
wie Koppelungskoeffizient und Koppelungsgrad verlieren für 
kontinuierliche Systeme ihre Bedeutung, ja direkt ihren Sinn. 
Will man streng vorgehen, so müssen jetzt die Gleichungen 
für die TM. als Randbedingungen für die Bewegungsgleichung 
der Luftsäule eingeführt werden. 

Die Untersuchung unserer Systeme ergab ein klares Bild 
von der Lage der Eigentöne, sowie von dem gesamten Me- 
chanismus. Jedes beliebige, aus Membranen und Lufträumen 
zusammengesetzte System läßt sich im wesentlichen auf unsere 
einfachen Modelle oder eine Kombination solcher Modelle 
zurückführen. 

Die weitere Untersuchung von Schwingungsystemen, die 
aus zwei oder mehreren Teilsystemen zusammengesetzt sind, 
und zwar die Untersuchung unter Abkehr von der rein punkt- 
förmigen Betrachtungsweise scheint lehrreich und lohnend 
zu sein. 


Für die Bereitstellung der in der vorstehenden Arbeit 
benutzten experimentellen Hilfsmittel haben wir der Helm- 
holtz-Gesellschaft und dem Elektrophysik-Ausschuß der Not- 
gemeinschaft der Deutschen Wissenschaft zu danken. Ebenso 
sind wir der Siemens & Halske A.-G. für Unterstützung der 
Arbeit zu Dank verpflichtet. 


(Eingegangen 9. September 1927) 
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4. Zur Theorie der Koppelschwingungen zwischen 


Telephonmembranen und Lufträumen;') 
von K. Schuster 


$1. Die Schwingungsgleichung der Telephonplatten und die 
Wellengleichung der Lufträume 

Im folgenden soll die strenge Theorie der in der vor- 
stehenden Arbeit von E. Waetzmann und K. Schuster 
„Koppelschwingungen kontinuierlicher Teilsysteme“ untersuchten, 
aus Telephonmembranen und Luftsäulen kombinierten Systeme 
gegeben werden. Von den dort behandelten Systemen I, II 
und III betrachten wir zunächst das System III, auf das I 
durch eine leichte Spezialisierung zurückgeführt werden kann. 


iL 
hy A 
6M 6 (2) 
T 7.M.(2) 
uf 
Fig. 1 


Das System III hat die in Fig. 1 skizzierte Gestalt. Das 
zylindrische Rohr X (Höhe A, Radius R) ist durch die beiden 
Telephonmembranen TM. (1) und TM.(2) begrenzt, die ihrer- 
seits durch die zylindrischen Gehäuse @ (Höhe h,) und G® 
(Höhe A,) vom Außenraum abgeschlossen sind. Die Wände 
des Rohres und der Gehäuse sollen als vollkommen starr be- 
trachtet werden. Die eine TM., etwa (1), möge unter der Ein- 
wirkung einer rein sinusförmigen Kraft stehen. Durch diese 
äußere Kraft wird das ganze System der beiden TM. und der 


1) Vorgetragen in der Sitzung des Gauvereins Thüringen-Sachsen- 
Schlesien der D. Phys. Ges. am 12. Juni 1927 in Breslau. 
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drei Lufträume zu erzwungenen Schwingungen erregt. Der 
Schwingungsvorgang bleibt infolge der Starrheit der Wände 
auf den von ihnen eingeschlossenen Raum beschränkt. Der 
Ansatz soll unter der Annahme der Zentralsymmetrie und 
ohne Berücksichtigung der Dämpfung gemacht, die TM. sollen 
als Platten aufgefaßt werden. 

Wir wollen von den Bewegungsgleichungen der beiden 
TM. ausgehen, wobei wir die von den Lufträumen herrühren- 
den Druckkräfte als äußere Kräfte betrachten. Nachher sollen 
die Vorgänge in den drei Lufträumen untersucht werden. Es 
wird sich zeigen, daß die Schwingungsgleichungen der TM. als 
besondere Randbedingungen für die der Wellengleichung unter- 
worfenen Vorgänge in den Lufträumen aufzufassen sind. 

Es sei y, die Elongation eines beliebigen Punktes der 
TM. (1), r sein Abstand vom Mittelpunkt der TM., e, die Dichte, 
D, die Dicke, Z, der Elastizitätsmodul, o, der Querkontraktions- 
koeffizient der TM. (1); P= P(r) gebe die Verteilung der er- 
zwingenden Kraft über die TM. an, » sei die aufgeprägte 
Frequenz, p der Druck im Rohr K, p'” der Druck im Ge- 
häuse G, p, der Normaldruck. Dann lautet die Schwingungs- 
gleichung der Telephonplatte (1): 


B 
+ Waren -AAy, = P-coswt 
— (P — Po) + (p™ — Pr) 


Dabei sind (p — p,) und (p'?— p,) die an den beiden Ober- 
flächen der TM. (1) herrschenden Uberdrucke, deren Verteilung 
und GréBe noch unbekannt ist. Da die TM. am Rande fest 
eingeklemmt ist, muß für r= R die Deformation und ihre 
Ableitung nach r verschwinden. Wir erhalten also die Rand- 
bedingungen: 


(2) r=R: y,=0; In 


(1) 


Ganz entsprechend gilt fir y,, die Elongation eines Punktes 
der TM. (2), die Differentialgleichung: 


wobei ¢,, D,, Z,, o, Dichte, Dicke, Elastizitätsmodul und Quer- 
kontraktionskoeffizient der TM. (2) und die doppelt überstrichenen 
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Ausdrücke die Überdrucke an ihren beiden Oberflächen be- 
deuten. Die Randbedingungen lauten analog zu Gleichung (2): 
r=R: 

Um die Beziehung zwischen den in den Schwingungs- 
gleichungen (1) und (3) vorkommenden Funktionen, den Elon- 
gationen der Platten und 
den Drucken, zu bestimmen, 
miissen die Schwingungen in r 
den Lufträumen untersucht 

Im mittleren Luftraum K 
sei s der Abstand eines be- 
liebigen Luftteilchens von der 
Platte (1), r sein Abstand von 
der Rohrachse (vgl. Fig. 2). 
Die Verschiebung q eines Luftteilchens kann aus einem 
Potential & abgeleitet werden: 


Zur Theorie der Koppelschwingungen usw. 


5:0 s-h 


Fig. 2 


Fir © = @(r, s, t) gilt die Wellengleichung 
(6) = a?- AO, 


wobei a die Schallgeschwindigkeit ist. Es müssen zunächst 
die Randbedingungen formuliert werden: Senkrecht zu einer 
starren Wand können die Luftteilchen keine Verschiebung er- 
leiden; daher müssen am Zylindermantel die Radialschwingungen 
verschwinden: 


(Ta) r=ak: 0. 


An den beiden Platten müssen die Längsschwingungen der 


Luftteilchen mit den Schwingungen der Plattenpunkte überein- 
stimmen: 


a@ 
(7b) s= 0: ds =Yı 
(Te) sak: oy, 


Diese Randbedingungen stellen die Verbindung zwischen den 
Schwingungen der Platten und der Lufträume her und können 
so als Ansatz der Koppelung betrachtet werden. 
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Ganz entsprechend sind die Verhältnisse in den Gehäusen, 
Im Luftraum G sei u der Abstand eines Luftteilchens von 
der starren Grundfläche des Gehäuses (vgl. Fig. 3). Dann 
existiert ein Verschiebungspotential ® = Di) (r, u, 2), das die 
Wellengleichung befolgt und den Randbedingungen genügt: 


(8a) r=k: 

(8b) u-0: 
8D 

(8c) u=h: 


Im Gehäuse G® sei z der Abstand eines Teilchens von der 
Platte (2) (vgl. Fig. 4). Für das der Wellengleichung unter- 


r=R 
r r 
— — — - > - — 
u 
u=0 u:h, Z=0 z=h, 
r=R r=R 
Fig. 3 Fig. 4 


worfene Verschiebungspotential — ®® (r, z, t) gelten die 
Randbedingungen: 


(9a) i =0, 

(9b) 

(9c) 20. 


Durch die Potentiale ® sind die Schwingungsvorgänge in 
den Luftriumen vollkommen bestimmt. An jeder Stelle ist 
die Verschiebung 


(10) q =— grad ® 
und die Druckschwankung 
(11) 


wo x das Verhältnis der spezifischen Wärmen der Luft und o 
die Dichte der Luft beim Normaldruck p, ist. 


] 
vorko 
(pP - 
poten 

(12a) 
(12b) 
| 
(12c) 
(12d) 
(12e) 
(12f) 
Setzi 

ein, 

(13) 

und 
Wi 
fass 
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§ 2. Die Formulierung des Randwertproblems 
Die in den Schwingungsgleichungen (1) und (3) der Platten 
vorkommenden unbekannten Funktionen y,, y,, (p — py), (P — Po) 


(p® — po) und (p® — p,) sind jetzt alle auf die Verschiebungs- 
potentiale zurückgeführt. Es ist 


(2b) y=— (52),.-- (=)... nach (7c) und (9b). 
Aus Gleichung (11) folgt: 


nach (7b) und (8c), 


(12c) (P — Po) = * Py) (4D), =0, 
(124) (p — Ps) = 2, 
(12e) (p'? — po) = # Py (4D), = n,» 
(p® 


Setzen wir die Ausdrücke in die Plattengleichungen (1) und (8) 
ein, so erhalten wir: 


(4.4.32), = Pecos oe 


rr as s= 


a@ 
— 8, D,- (5° xp, (4D), =» 


Dabei ist 
v=1,2 


und 4 das Symbol des Laplaceschen Operators für die 
Koordinate r: 
6? 1 @ 
tr 

Die Plattengleichungen haben also die Gestalt besonderer 
Randbedingungen für die Verschiebungspotentiale angenommen. 
Wir können jetzt das ganze Problem folgendermaßen zusammen- 
fassen: Es sind die Funktionen ®, ®® und ®® zu bestimmen. 
Sie sollen der Wellengleichung genügen: 


529 
— xp, (4D), <9. 
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| b 
DV = a? 4D», 
| = a2 4 
und die Randbestimmungen erfüllen: 


(5a) s=0, u=h: 


(14) 


nach (12a), und 


a@ 
14,4, 5, = P-cosmt—xp,A® 


+«p,4®” nach (13), 


(15c) s=h, z=0 se = en nach (12b), und 
ö® 

nach (13), 

Da 

(15e) r=R: “a = — = 2s = 0 nach (7a), (8a), (9a), 

(15f) «=0: =0 nach (8b) 

(lög) z=h,: =0 nach (90), 

(15h) s=0, 22 =0 nach 

(Ibi) s=hr=R: nach (4), 


Da die Dämpfung vernachlässigt ist, finden die Schwin- 
gungen einer bestimmten Frequenz in allen Teilen des Systems 
in gleicher Phase statt; wenn-es sich dabei um erzwungene 
Schwingungen handelt, haben sie die gleiche Phase wie die 
äußere Kraft. Wir werden daher setzten: 


® = wfr,s)-coswt, 
(16) DV = y(r,u)-coswt, 
D® = z)-cos@t. 


Es sollen zunächst die freien (P(r)=0) und dann erst die er- 
zwungenen Schwingungen behandelt werden. Im ersten Fall 
ist die Frequenz w eine noch zu bestimmende Größe, nämlich 
ein Eigenton des Systems; im zweiten Fal) ist m die gegebene 
aufgeprägte Frequenz. 


(18) 


gleich 
übers 
aus & 
es Vi 
blem: 
Ritz 
werd 
Polyı 
(17) 
Die 
80 W 

und 
Koef 

daB 

gut“ 

Glei 

(19) 
so d 

(20) 
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Für die Raumfunktionen y, yw, ıy'® gelten dann die 
gleichen Beziehungen wie fir ©, ©, @®, wenn man für ® 


überall (— ©?) einsetzt, da sich der zeitliche Faktor cos wt 
aus allen Gleichungen heraushebt. 


§ 3. Die Methode von Ritz 


Wegen der Kompliziertheit der Randbedingungen erscheint 
es von vornherein aussichtslos, die Eigenfunktionen des Pro- 
blems in exakter Form zu finden. Es soll daher die von 
Ritz!) angegebene Näherungsmethode zur Integration benutzt 
werden. Wir setzen für die gesuchten Funktionen y, y?, y® 
Polynome von besonderen Funktionen an: 


0,n 
0, n 
(17) = Sa,w,®, 
0,n 
= >, v,, 


Die Funktionen »,(r,s\, w,(r,s), w,2(r,u), v,®(r,z) wollen wir 
so wählen, daß die Funktionen ®, OD”, D® die Wellengleichung 
und möglichst viele Randbedingungen streng erfüllen. Die 
Koeffizienten a,, b, werden wir dann so zu bestimmen suchen, 
daß die noch übrig bleibenden Randbedingungen „möglichst 
gut“ erfüllt werden. Wir setzen: 


cos a, (s — h) : 


COS a, 8 
(18) „tr 
_ (% — hy) -J, (8 r). 


sin a, hy 


Dabei seien die £,(v = 0,1,...n) die n + 1 ersten Wurzeln der 
Gleichung 


(19) 4, (8, = 0, 
so daß im besonderen f, = 0 ist. Weiterhin sei 
(20) k=" ud =k 


1) W. Ritz, Ann. d. Phys. 28. 8. 737. 1909. 
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so daß «, =k ist. Ferner wollen wir für die a, und die J, 
die Beziehungen vorschreiben: 


0,n 
(8, 2) 0, 
(21) 


0,n 
(8, = 0. 


Auf diesen Ansatz (17) bis (21) kommt man, indem man die 
Elementarlösungen der Wellengleichung in Zylinderkoordinaten 
aufsucht und dabei die an den starren Wänden geltenden 
Randbedingungen berücksichtigt. Die Gleichungen (21) stellen 
die Bedingung dafür dar, daß die Ränder der Telephonplatten 
in Ruhe bleiben. Bei der angegebenen Wahl der Funktionen 
w, erfüllen die ®, die Wellengleichung (14) 
und die Randbedingungen (15) außer den beiden (15b) und (154), 
die den Schwingungsgleichungen der Platten entsprechen. 

Es handelt sich jetzt noch darum, das System der a,, }, 
so zu bestimmen, daß diese beiden Randbedingungen „mög- 
lichst gut“ erfüllt werden. Das geschieht nach Ritz, indem 
wir mit unserem Näherungsansatz in ein Variationsprinzip ein- 
gehen. Wir wollen aus den Plattengleichungen (P= 0): 
(22) & D, e,44y, = —(p—Pp) (rdr d6)-dy, 

& + = +(P—Po) —(P—Py)| (rar dy, 
die allein durch unseren Ansatz nicht befriedigt werden, das 
d’ Alembertsche Prinzip herstellen, indem wir mit den Flächen- 
elementen (rdr d6) der Platten und mit virtuellen Verschiebungen 
öy, und öy,, die beliebige Funktionen von r sind, erweitern, über 
die Platten integrieren und die beiden Gleichungen addieren: 

R R 


[in [ Ady, -dy, rar 
6 


R R 
+ fo = pa) dy, 7 dr — (p® = p,)- dy, rar 
(23) 0 0 


R R 
D, dy, 7dr + [449 rar 
0 0 


R R 
| - [ep -öy,rar + [(p? =p,)-dy,rdr =0. 
0 0 


(2 


aus 
dri 
Pri 
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Wi 
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Setzen wir jetzt für die Funktionen y,, 7,,(p—p,) usw. die 
aus dem Polynomansatz für yw, w'’, w® sich ergebenden Aus- 
drücke ein, so stellt die Gleichung das d’Alembertsche 
Prinzip nicht nur für die Plattenschwingungen, sondern für 
den Vorgang im gesamten Schwingungssytem dar, weil die 
Wellengleichung und die übrigen Randbedigungen durch den 
Näherungsansatz streng erfüllt werden. Zu variieren ist jetzt 
nach den Koeffizienten a,, 6,, über die allein ja noch nicht 
verfügt ist. Das vorher noch unbestimmte „möglichst gut“ ist 
jetzt dadurch genau definiert, daß der Vorgang nach dem 
d’Alembertschen Prinzip verläuft, wobei allerdings die Willkür 
der virtuellen Verschiebung durch den Näherungsansatz ein- 
geschränkt ist. 
Aus den Gleichungen (12) und (16) bis (18) ergibt sich: 


Ow). 


=— cosmt- (=)... 
(24) | 
= cosm@t- > %%(8,7) 


v 


0,2 
dy, = coswt- > 


Auf gleiche Weise kénnen wir alle in (23) vorkommenden 
Funktionen darstellen. Setzen wir die gefundenen Werte in 
(23) ein, so nimmt diese Gleichung nach Vertauschung der 
Summationen mit der Integration und unter Berücksichtigung 
der Orthogonalitätsbeziehung 


R 
fra für 
und 
2 
frard2(,r) = 


0 


folgende Gestalt an: 


b, 
die 
ten 
den 
llen 
sten 
nen 
(14 
5d), 
lem 
ein- 
Yı 
Yo = 
1en- 
gen 
iber 
ren: 
). u 
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0? + 0, B+ 00° 


a, a, sin a, h 


(25) 


Diese rn, wollen wir u schreiben: 
(26) A,-da, + - B,+ db, =0. 


Da nun de Variationen da,, 3b, eine virtuelle Verschiebung 
darstellen sollen, miissen sie a den Randbedingungen in Kin- 
klang stehen. Die den Randbedingungen (15h) und (15i) ent- 
sprechenden Relationen (21) müssen also auch bei der Variation 
erhalten bleiben. Es muß also sein: 


0,n 
(27) 


> 5b, = 0 u 


Die Gleichungen (26), (27) stellen eine Extremumsaufgabe mit 
Nebenbedingungen dar. Nach dem Lagrangeschen Multipli- 
katorenverfahren ergibt sich dann: 
4,+i=0 

| B,+u=0 
Wir setzen jetzt zur Abkürzung: 
(29) | 4, J, (8, R)-(a,)- a, +4, (8, 2)-(k,) b 
B, = + R)-(0,)-b, 
Es ist also nach (25) und (26): 

(a,) &D, eo” + be | ß,* + 0 0°. + eotg a, ) 


ay 


(30) (,) = — + 48,4 + + 


ay a, 
1 


y=0,1,...n. 


(k,) =— - 


Du 


(1) 
Fil 
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Durch Einsetzen in (28) erhalten wir: 
(8, B)+(h,)-a, + (8, R)-(b,)-b, + 0 


Fügen wir noch die beiden Gleichungen (21) hinzu. 


0,n 
HGR) =0, 


0, n 
=0, 


so haben wir für die 2n +4 Unbekannten a,, b,, A, u ein System 
von 2x + 4 linearen homogenen Gleichungen erhalten. Wenn 
ng § dieses Gleichungssystem lösbar sein soll, muß die Determinante 
n- | der Koeffizienten verschwinden. Da diese nach (30) Funktionen 
it | der Frequenz w sind, erhalten wir eine Gleichung von der Form 


on (32) F(o) = 0. 
Nach leichter Umformung der Determinante nimmt diese 
Gleichung die Gestalt an: 
(a,)+ (a,) («,)}— (a,) (hg) + (A,) (k)— = (k,) 
(a,) (a). (A,) 
| \ 
(hy) (Ay) (b,) (,) 
Na) 


An den freien Stellen stehen Nullen. Die Wurzeln der Glei- 
chung /(w) = 0 liefern eine Reihe von diskreten Frequenzen, 
die Eigentöne des Schwingungssytems. Diese Eigenfrequenzen 
werden mit um so größerer Annäherung und Vollzähligkeit er- 
halten, je größer die Zahl n gewählt wird. Jedem Eigenton 
entspricht ein bis auf einen gemeinsamen Faktor bestimmtes 
System der a,, b,. Infolgedessen können auch die Eigenfunk- 
tionen w, y'?, w® mit beliebiger Annäherung berechnet werden. 

Es sollen jetzt kurz die erzwungenen Schwingungen unter- 
sucht werden. Es ist also P= P(r) eine gegebene Funktion, 


a} 
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die die Verteilung der erzwingenden Kraft längs der Platte (1) 
angibt, und w die gegebene aufgeprägte Frequenz. Wir ent- 
wickeln P(r) in der Form: 


0,n 
P(r) > %%(8,7)- 


Das ist natürlich nur mit einer gewissen Annäherung möglich, 
da n eine endliche Zahl ist. Dann können wir in der gleichen 
Weise verfahren wie bei den freien Schwingungen. Jetzt er- 
halten wir aber für die Größen a,, b, ein inhomogenes lineares 
Gleichungssystem. Dieses ist nur lösbar, wenn die Determinante 
der Koeffizienten von Null verschieden ist. Die Determinante 
ist die gleiche wie bei den freien Schwingungen. Wir er- 
halten also: F(o) + 0. 


Es werden also die Stellen = 0 von der Lösung aus 
geschlossen, d. h. es darf die aufgeprägte Frequenz w nicht 
mit einer der Eigenfrequenzen des Systems zusammenfallen. 
Das ist ganz natürlich, da es sich um ungedämpfte erzwungene 
Schwingungen handelt, die an den Resonanzstellen F(a) = 0 
unendlich große Amplituden annehmen würden. In jedem 
Fail, wo F(o) +0 ist, können die a,, 6, durch die c, aus- 
gedrückt werden, wodurch der Schwingungsvorgang vollkommen 
bestimmt ist. 


$ 4. Die Bestimmung der Eigentöne 
Im folgenden sollen die gewonnenen Resultate in der 


ersten Näherung der Ritzschen Methode untersucht werden. 
Die Frequenzengleichung (33) lautet für n=1: 


(34) ()+(a) (A) + (A) 

(ho) + (Ay) |- 
Wir wollen jetzt die vereinfachende Annahme machen, dab 
die beiden Platten und die beiden Gehäuse völlig gleichartig 


sind. Dann wird 
D, =8,D,=eD, 
= Cc, 
hy = hy. 
Infolgedessen ist nach (30) (a,) = (6,) und wir erhalten als 
Frequenzengleichung: 


(35) = + + (A,)] 
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oder, wenn wir die Werte aus (30) einsetzen: 
— 2:Dw? + ow? 


„[eotgk h cotg a, h cotg kh, cotg a, h, ) 
( k + ay k + a, 


(36) 


1 1 
+ 
Wenn der Radius R klein gegen alle vorkommenden Wellen- 
längen ist und die Höhen A und %, nicht sehr klein gegen R 
sind, können die Glieder mit «, v vernachlässigt werden. Weiter- 
hin seien die Höhen A, =A, der Gehäuse klein gegen die 
Wellenlänge. Dann wird 


2, cotgkh, _ Qe 
(87) 09 k h, 


und wir erhalten: 

A kh 1 
Setzen wir noch zur nae 


so wird endlich 


Dabei ist offenbar w = » der Eigenton der mit dem Gehäuse 
verbundenen Platte, wenn von der Seite des mittleren Luft- 
raumes her keine Kräfte auf die Platte wirken. Da die Un- 
bekannte © auch im Argument des cotg und sin vorkommt 
(denn k ist gleich w/a), läßt sich die Gleichung (40) nicht 
nach » auflösen. Um den Zusammenhang zwischen der Lage 
der Eigenténe w und der Höhe A zu übersehen, empfiehlt sich 
eine graphische Darstellung. ') 

Wir wollen jetzt die Frequenzengleichung für das 
Schwingungssystem I aufstellen; dieses geht aus dem bisher 
behandelten Fall III dadurch hervor, daß die eine Telephon- 
platte, etwa (2), durch eine starre Wand ersetzt wird. Dann 
werden in unseren Gleichungen alle 5, gleich Null und als 
Frequenzengleichung ergibt sich in der n-ten Näherung: 


1) E. Waetzmann u. K. Schuster, a. a. O. Fig. 13. 
Annalen der Physik. IV. Folge. 84. 85 
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(a) +(a) —(a)— —(@,) 


(41) (4) 
(a,) (a,) | 

Für n=1 wird 

(42) (a,) + (a,) = 0. 


Setzen wir in diese Gleichung die Werte aus (30) ein und 
treffen die gleichen Vereinfachungen wie oben beim System III, 
so erhalten wir: 


Hieraus sieht man sofort, daß bei geeigneter Wahl der Hohe h 
ein Eigenton w des Systems mit dem Eigenton r der mit ihrem. 
Gehäuse verbundenen Platte zusammenfällt. Wenn w=n 
werden soll, muß der cotg verschwinden, also sein Argument 
gleich einem ungeraden Vielfachen von 2/2 sein. Das ergibt 
im einfachsten Fall: 


a a 2 


Infolgedessen wird h = +, wenn A, die dem Eigenton n ent- 


sprechende Wellenlänge ist. Auf diesen Spezialfall werden 
wir im nächsten Paragraphen noch einmal zurückkommen. 


$5. Die Amplitudenverteilung 


Es soll jetzt kurz auf die einem Eigenton zugeordnete 
Amplitudenverteilung längs der Telephonplatten und in der Achse 
des Luftzylinders eingegangen werden. Sie ist im ersten Fall 


gegeben durch (37), baw: (32), im zweiten durch 
5), Kr In der ersten Näherung der Ritzschen Methode ist 


cosk(s — h) cos ks cos a, (s — h) 
(44) | ~ S 


(A, r) + b, I, (A; r). 


a, sina, h 


Infolgedessen ist 
ö 
| a) + a, + J, (8, 7) 


(45) 
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Da nach (21) die Beziehungen gelten: 


146) a, +a, +J,(8, R)= 0 
by + b, = 0, 
erhalten wir als Profile der beiden Telephonplatten: 


| 


Da in der ersten Näherung das Profil durch das Verhältnis a,/a, 
bzw. b,/d, vollkommen bestimmt ist, diese Quotienten aber 
durch Gleichung (46) festgelegt sind, entspricht der ersten 
Näherung die Annahme, daß die Telephonplatten unabhängig 
von allen äußeren Verhältnissen ein bestimmtes, für alle Eigen- 
tine gleiches Profil haben. 


In der Achse des zylindrischen Rohres wird 


(47) 


cos k(s — h) ma by cosks 
148) Yr=0=— ksin kh a, ksinkh } 
cos a, (s — h) 6, cosa,s 
1 a, sin a, h @,sina,h 


Die Verhältnisse 4,/a, und 4, /a, ergeben sich aus dem Gleichungs- 
system (31): 
149) bo __W+@), 

% (ko) + (le) 


Nun gilt bei gleicher Beschaffenheit der beiden Platten für 
einen Eigenton nach (35): 


(a,) + (a,) = F [(#,) + (4,)]- 


Folglich ist 
b b 

(50) = =+1. 

Wir wollen je nach dem Vorzeichen in Gleichung (35) zwei 
Gruppen von Eigentönen unterscheiden, ®, (wenn in (35) 
das — gilt) und w_ (wenn in (35) das + gilt). Dann lehrt 
Gleichung (50): Bei jedem Eigenton des Systems III schwingen 
die beiden, als völlig gleichartig angenommenen Telephonplatten 
mit gleicher Amplitude und zwar bei einem Eigenton w, in 
gleichem, bei einem Eigenton w_ in entgegengesetztem Sinn. 

35* 
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Infolgedessen wird 


—_, . [cosk(s h)Feosks _ 1 
| kein kh B) 
a, sin a, h 


und wir erhalten als Amplitudenverteilung in der Rohrachse 
bei einem Eigenton w,: 


a cos-* (28-1) 1 cos (25-1) 
cos ow cos th 
bei einem Eigenton w_: 
sinh sinh 


Man sieht hieraus, daß im ersten Fall in der Mitte des Luft- 
zylinders ein Schwingungsbauch, im zweiten Fall ein Schwingungs- 


knoten liegt. 


S=0 S=h4 
Fig. 5 
Ganz entsprechend erhalten wir bei System I, wo die 
Platte (2) durch eine starre Wand ersetzt ist, folgende Ampli- 
tudenverteilung in der Rohrachse: 


Oy sink(—h) 1 sin a,(s — h)\ 
(64) (2 a, { sinkh Jo (8, R) sina,h ) 
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An der starren Wand (bei s = 4) liegt naturgemäß immer ein 
Schwingungsknoten. 

Als Beispiel diene der am Ende des vorigen Paragraphen 
erwähnte Spezialfall: Wenn beim System I die Höhe A des 
Zylinders gleich der Viertelwellenlänge des Eigentons n der 
mit ihrem Gehäuse verbundenen Platte ist, wird ein Eigen- 
ton ® des Systems gleich n. Für diesen wollen wir die 


Amplitudenverteilung untersuchen. Da jetzt hh = + ist, er- 


gibt sich aus Gleichung (54): 
1 sin a, (s —h 

(55) (3%), {sin} (+ 1) + ze h 
Wenn der Radius # klein gegen die Wellenlänge ist, kann 
man sich leicht davon überzeugen, daß der zweite Summand 
nur in unmittelbarer Nähe der TM. einen merklichen Beitrag 
zur Amplitudenverteilung liefert. Dieser Fall ist in Figur 5 
dargestellt. 


$ 6. Behandlung des Problems im Anschluß 
an die Punktmechanik 

In § 4 wurde die Bestimmung der Eigentöne in der ersten 
Näherung der Ritzschen Methode und unter weiteren ver- 
einfachenden Annahmen beschrieben. In § 5 wurde gezeigt, 
daß hierbei den Telephonplatten bei allen Eigentönen das gleiche 
Profil zukommt. Das besagt natürlich, daß die erste Näherung 
nur für die Umgebung des Grundtones ausreicht. Es liegt 
nahe, unter Beibehaltung der Beschränkung auf kleine Frequenz- 
bereiche die Frequenzengleichungen auf einem einfacheren 
Wege als nach Ritz abzuleiten, indem man die Vereinfachungen 
schon in den Ansatz einführt. Dann verzichtet man natürlich 
von vornherein auf eine genaue Berechnungsmöglichkeit, wie 
sie die Ritzsche Methode gibt. Aber in manchen Fällen, 
80 bei dem bisher noch nicht behandelten System II, führt 
das abgekürzte Verfahren leichter zum Ziel. 

Wir wollen jetzt von vornherein das Profil als eine Eigen- 
schaft der TM. auffassen und als bekannt ansehen. Das Profil 
der TM. (1) ist durch das Verhältnis der Koeffizienten a, zu- 
einander gegeben. Das System der a, bestimmte nicht nur 
die Form, sondern auch die Größe der Durchbiegung der 
TM. (1). Da jetzt aber nur die Form, das „Profil“, gegeben 


_ 
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sein soll, wollen wir die Verhältniszahlen der a, einführen; 
trotz der abweichenden Bedeutung soll für diese neuen 
Koeffizienten die Bezeichnung a, beibehalten werden. Die 
Verhältniszahlen a, seien so gewählt, daß 


0,0 


(56) Da, = 1 e 


ist. Ganz entsprechend wird das Profil der TM. (2) durch die 
Verhältniszahlen 5, festgelegt, fir die die Gleichung gilt: 


(57) 


Wir untersuchen zunächst das System III, das aus zwei durch 
ein zylindrisches Rohr verbundenen Telephonen bestand. 
Sind , und z, Zeitfunktionen von der Form: 


(58) zr, = Acosot, 
x, = Bcoswt, 
so können wir nach (16), (17), (18) das Verschiebungspotential 
im mittleren Luftraum folgendermaßen schreiben: 
0,0 
inch" 20 8,7) 
(59) 


8a, 


Dabei sind die Größen A-a, und B- d, die früheren Koeffi- 
zienten a, und 6,. Es läßt sich nun sofort zeigen, daß z, 
und z, die Mittelpunktselongationen der TM. darstellen. Die 
Elongationen beliebiger Membranpunkte sind nach (12)und (59): 


0,0 
= 
(60) 


Die Funktionen f,(r) und f,(r) sind dabei die Profile der 
TM. Die Mittelpunktselongationen erhalten wir, wenn wif 
r = 0 setzen: 


(61) 


(6 


= 
| 
spr 
62 
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0,0 
x, nach (56) 
(61) 


0, @ 
(Ya) =0 =. =2,, nach (57) 


x, und z, sind also wirklich die Mittelpunktselongationen der TM. 
Die Potentiale in den Gehäuselufträumen lauten ent- 
sprechend: on 


cos a, (x — hy) 
ay (% — 


Da die a,, 5, gegeben sind, hängt der Zustand des ganzen 
Systems nur noch von den Größen z, und z, ab. Wir werden 
sie daher als Koordinaten benutzen. Wir wollen jetzt auf das 
System der beiden TM. ohne die Lufträume die Lagrangeschen 
Gleichungen anwenden und die von den Lufträumen herrühren- 
den Druckkräfte als äußere Kräfte ansehen: 


d (oT oT 

(64) 
dt \ da, 02, 


Dabei sind 7 und U kinetische und potentielle Energie 
der beiden TM. Die Kräfte X, und X, sind durch die bei 
einer virtuellen Verschiebung geleistete Arbeit ö’ A definiert: 


(65) 
Die kinetische Energie der beiden TM. ergibt sich durch 
Integration über die Energien der einzelnen Membranpunkte zu 


R R 
0 0 
R 
1 
= 54": ar) 
0 


R 
0 
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Da die Profile f,(r) und f,(r) gegeben sind, sind auch 
die in den runden Klammern stehenden Ausdrücke, die man 
als Schwingungsmassen der TM. bezeichnet, bekannt; wir 
nennen sie m, und m,. Dann wird 
(66) T = }m, &?-+ m, #,). 

Die beiden TM. können also rein formal als Massen- 
punkte betrachtet werden. Entsprechend sei die potentielle 
Energie der TM. 

(67) 
wobei c, und c, die Elastizitäten der TM. sind. Durch Ein- 
setzen in (64) erhalten wir als Schwingungsgleichungen: 
(68) 

m, +c, 2, = X,. 
Es handelt sich jetzt noch darum, die auf die TM. wirkenden 
änßeren Kräfte X, und X, zu bestimmen. 


$7. Die von den Lufträumen herrührenden Druckkräfte 

An äußeren Kräften wirken auf die TM. (1) die er- 
zwingende Kraft P(r)-cosw¢, vom mittleren Luftraum der 
Druck (p — p,), vom Gehäuse der Druck (p — p,); auf die 
TM. (2) wirken von beiden Seiten her die Drucke (p — p,) 


und (p?—p,). Um aus diesen Drucken die Kräfte X, und 
X, zu erhalten, müssen wir die bei einer virtuellen Verschie- 
bung der beiden TM. geleistete Arbeit 0’ 4 bilden: 

R 


2ncoswt. {rdr P(r). Oy, 
0 


R 

(69) | + rdr(p®—p,)-dy, 
0 

+2n “fr dr(p 
0 


R 
— 2m- [rdr(p® — py)- 
0 


(72) 


(73 


Int 


e 


d 
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Da die Profile bekannt sind, werden die Variationen nach 
(60) jetzt einfach 


(0) 


so daB wir erhalten: 


f 
0’ A = 62, [om f [rar Po) 
0 


=f, 
Dy, = 


—2a-frdr(p—p) 
0 


R 
m) | dr(p® — py) 
R 


R 
—2n- frdr(p® — p,)- val 
0 
| = X62,4+ X,62,. 
Durch Vergleich der Faktoren von dz, und dz, erhalten 
wir die Krafte X, und X, 
Nun ist nach (12) und (59): 


(p — = a (4D)-0 = 0 (@),=0 


0, @ 


0,0 
1 
+ b, %0(8,7) . 


In gleicher Weise lassen sich auch die anderen Drucke 
darstellen. Außerdem waren die Profile 


fiir) = Sa, 


0,0 
f(r) = 2b, J, (8,7) 


Setzen wir diese Werte in (71) ein, so treten unter den 
Integralzeichen überall Doppelsummen auf. Infolge der 


(13) 
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Orthogonalität der J,(8,r) bleiben aber nur die Integrale über 
die rein quadratischen Glieder stehen, so daß wir erhalten: 


(A= Oz, 
0,0 
0,» 
1 
0, a 
(74) v 
0, 
1 
+dz,° I- a,-sina,h J," (8, B) 
0,00 
h 
0,0 
| + = b2. enh 
Dabei ist gesetzt 2 


Q= 2a [rar 


Die geschweiften Klammern stellen die Kräfte X, und X, 
dar. Wir setzen zur Abkürzung: 


X, = Qcoswt+ K, —K,,+G, 

= K,, + K, + G,. 

Dann ist Qcoswt die erzwingende Kraft; (K, — K,,) und 
(X, — K,,) stellen die vom mittleren Luftraum auf die beiden 
TM. ausgeübten Kräfte dar; G, und G, rühren von den Ge- 
häusen her. Unsere Schwingungsgleichungen (68) lauten dann: 
(76) + K,, =Qcoswt, 

m, #, + ¢,2,—G,+K,,—K, =0. 

Wir wollen jetzt die von den Lufträumen herrührenden 
Kräfte noch weitgehend vereinfachen, indem wir von den sie 
darstellenden Summen nur das erste Glied berücksichtigen. 
Die übrigen Glieder können — zwar nicht gegen das erste, 
denn dieses kann bei geeigneter Höhe Ah gleich Null werden, 


(75) 


ver 


wohl 
nach 
kom 
auße 
nehı 
die 
Dan 
forn 
(17) 
wir 
sie 
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Da 
(79 
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wohl aber gegenüber den Trägheitskräften der TM. — ver- 
nachlässigt werden, wenn der Radius & klein gegen alle vor- 
kommenden Wellenlängen und die Höhen A, h,, A, nicht 
außerordentlich klein gegen R sind. Ferner wollen wir an- 
nehmen, daß die Höhen A, und A, der Gehäuse klein gegen 
die Wellenlänge sind, daß also Ah, <1 und kh, <1 ist. 
Dann lassen sich die Gehäusekräfte in folgender Weise um- 
formen: 


G, =i, o’n R* ga,? 


1 Ph, 
n R* 9 a* 0,2 
(77) 
G,=—2,° b, 


Da die Faktoren von x, und z, frequenzunabhängig sind, 
wirken die Gehäuse als reine Zusatzelastizitäten; wir wollen 
sie daher mit c, bzw. c, zusammenfassen, so daß c, und c, 
von jetzt ab die Gesamtelastizitäten der mit ihren Gehäusen 
verbundenen TM. darstellen. 


Ferner wollen wir zur Abkürzung setzen: 
(18) 


1 
@Ma- n R* 0 a, by . = 
Dann wird 
(79) =h, +2; —Ky, = 
— Ky, = 23 —K, =h -2,. 

Die Schwingungsgleichungen (76) nehmen dann die Ge- 
stalt an: 

(80) +(e, +) 2, — = Jeoswt 
m, + (cg th), — = 0. 

Um die entsprechende Gleichung fir das SystemI zu 
erhalten, brauchen wir nur, da hier die TM. (2) durch eine 
starre Wand ersetzt ist, z, = 0 zu setzen. Wir erhalten also: 
(81) m, #, + (c, + k,)2, = Qcoswt 

Die Gleichungen (80) und (81) unterscheiden sich von den 
gewöhnlichen Massenpunktsgleichungen nur dadurch, daß die 
Koeffizienten A,, k,, k,, frequenzabhängig sind. 
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§ 8. Die Bestimmung der Eigentöne 

Für freie Schwingungen gehen die Gleichungen (80) durch 
den Ansatz (58) über in: 
(82) (-o’m + ¢, +h) +2, — =0 

+(—w?m, +, + h,)- 27, = 0. 

Soll das Gleichungssystem lösbar sein, so muß die Deter- 
minante verschwinden; das liefert die Gleichung für die Eigen- 
tine. Wir wollen wieder gleiche Beschaffenheit der beiden 
Telephone annehmen. Dann wird: 


m=m=m, 
kh=h, 

und wir erhalten als Frequenzengleichung: 

(83) —w?m+ce+h = 


Setzen wir die Werte aus (78) ein, so wird 
e, nR’oa-a, 1 
(cotghh 


m 
Nennen wir das Verhältnis von Schwingungsmasse zu 
Gesamtmasse der TM. u, so ist m=p-ah?eD. Infolge- 
dessen wird 
m D u eD 
Ferner sei — = n?. Dann ist n wieder der Eigenton der mit 


ihrem Gehäuse verbundenen TM. und wir erhalten: 


u sD 

Gleichung (85) unterscheidet sich von (40) nur dadurch, 
daß der dort auftretende Faktor } jetzt durch die — nur vom 
Profil abhängige — Größe a,?/u ersetzt ist. Das liegt natür- 
lich daran, daß vorhin das Profil sich zwangsläufig aus der 
Rechnung ergab, während es jetzt willkürlich festgesetzt ist. 
a, hat dabei die Bedeutung des Durchbiegungsvolumens der 
TM. im Verhältnis zu dem Volumen des Zylinders, der den 
gleichen Radius wie die TM. und ihre Mittelpunktsamplitude 
zur Höhe hat. 

Für das System I ergibt sich aus (81) die Frequenzen- 
gleichung: 
(86) —a’m+e+h, =0 


nu 
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oder nach Einsetzen der Werte: 

(87) — nt cotg hh, 

in Analogie zu Gleichung (43). 


§ 9. Das Schwingungssystem II 
Wir wollen jetzt das in Fig. 6 skizzierte System II unter- 
suchen. Der Luftraum G®” ist jetzt ein enges Rohr vom 
Radius #,, dessen Höhe A nicht mehr klein gegen die Wellen- 
länge sein soll. Zwischen den Lufträumen X und G® befindet 
sich jetzt keine TM. mehr. Wir können uns das System in 
folgender Weise aus dem System III entstanden denken: Das 


IL 


h, h c 
@) N (2) 
CG x 6 
ZM\K | 
Fig. 6 


Profil f,(r) der TM. (2) sei in der Weise verändert, daß es im 
Bereich von r=, bis r=R gleich Null ist, so daß in 
diesem Bereich die TM.(2) jetzt eine starre Wand darstellt. 
Über die Schwingungsform im Bereich von r = 0 bis r = R, 
wollen wir keine Annahmen machen. Da sich nun beim 
System II zwischen den Räumen X und G® gar keine TM. 
befinden soll, müssen wir die TM.(2), mit der wir bisher 
operiert haben, als masse- und elastizitätslos annehmen, also 
m =0 und ,=0 setzen. Die Funktion f,(r) stellt also in 
Wirklichkeit das Profil der in der Trennungsebene von X und 
G® schwingenden Luftteilchen dar. 

Dem neuen Radius A, entsprechend lautet jetzt das Ver- 
schiebungspotential ®® im Raume G® analog zu (63): 

On 


„( -h 
(88) 


Dabei ist z, die Mittelpunktselongation der in der Tren- 
nungsebene von K und G® (also s=h, z= 0) schwingenden 
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Luftplatte. Die 5, werden durch die Beziehung J, (6, R,) = 0 
bestimmt, die 7, ergeben sich aus den Gleichungen y,? = k? —ö? 
(vgl. (19) und (20)). Hieraus folgt für das Profil der Luft- 
platte im Bereich von r=0bisr=AR,: 


0,n 
(89) f(r) = Se, J, (4,7); 
im Bereich von r=0 bis r = KR wird f,(r) wie früher dar- 


gestellt: 0,n 
(90) = = b, J, (8,7). 


Aus dem Potential ®® berechnet sich der Druck auf die Luft- 
platte vom Raume G® her: 


(91) (p® — po — Po) = 0 (D?), =0° 
Wir erhalten daher als Kraft, die vom Luftraum G® 
herrührt: 
RB, 
6, = 2n-[rarip® — Po) fe (r) 
0 


0,0 
un Hy +m .J,2(6, 
(92) | d 
= cotz kh)-z. 
2 


| + 72 
Die Größe y,, die die höheren Glieder der Summe darstellt, ist 
frequenzunabhängig, wenn A, klein gegen die Wellenlänge ist. 

Für das System II gehen nun nach den obigen Betrachtungen 
die Schwingungsgleichungen (76) über in die Gleichungen: 
(93) | m, #,+¢,2,—G,—K,+ K, = Qcosat, 

— G+ — K, =. 

Mit Ausnahme von G,, was soeben neu hergeleitet wurde, 
haben alle Buchstaben die gleiche Bedeutung wie beim 
System III. Nur wollen wir jetzt in der zweiten Gleichung bei 
den Kräften G, und K, die höheren Glieder der Summen mit 
berücksichtigen, weil hier keine Trägheitskräfte vorhanden 
sind, denen gegenüber die entsprechenden Größen in der ersten 
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Gleichung vernachlässigt werden können. Die vom Gehäuse G™ 
herrührende Kraft G, vereinigen wir wie früher mit dem 
Glied c,z,. Die übrigen Kräfte lauten dann (vgl. (79)): 
(94) -Kehz; — Ky, 
— Ky = - = 9,2, + 
Die Größe x, %, stellt dabei die höheren Summanden von — K, 
dar, x, ist frequenzunabhangig, wenn R klein gegen die Wellen- 
lange ist (vgl. (92)). Unsere Schwingungsgleichungen (93) 
nehmen dann die Gestalt an: 
(95) m, #, — hz, = Qcosamt, 
(hy + + + 79) — = 0. 
Für freie Schwingungen erhalten wir daher: 
(96) (— m, + + — 2, = 0; 
+ (hy + 9, — [x = 0. 

Das Verschwinden der Determinante liefert die Frequenzen- 

gleichung: 

Wir wollen jetzt annehmen, daß die Höhe A, des mittleren 
Luftraumes klein gegen jede vorkommende Wellenlänge ist. 
Dann vereinfachen sich die Gleichungen (68) — dort ist die 
Höhe des mittleren Luftraumes h genannt — in folgender Weise: 

, 
98) 


a R* oa°-a,b, : 
12 


Infolgedessen ist A?, = A, A, und wir erhalten: 


m 
Setzen wir für g, und A, die Werte aus (92) und (98) ein, 
so wird: 


R?b,?-a 
(100) ctgkh= — . 


2 


mM; 
+ 
m; 


552 K.Schuster. Zur Theorie der Koppelschwingungen usw. 


Da 4, und c, das Durchbiegungsvolumen der Luftplatte dar- 
stellen, und zwar 4, im Verhältnis zum Zylinder vom Radius A, 
c, im Verhältnis zum Zylinder vom Radius R,, wobei die Höhe 
beider Zylinder gleich der Mittelpunktsamplitude der Luftplatte 
ist, so ergibt sich: 


b _ R,* 
(101) 
Setzen wir noch zur Abkürzung 
+7; 


so erhalten wir als Frequenzengleichung für das System II: 


(102) cotg kh Bho 


Zusammenfassung 


An einigen speziellen Systemen, die sich zur theoretischen 
Behandlung eignen, werden die Schwingungen gekoppelter 
Kontinua untersucht, nämlich an einigen Kombinationen von 
Telephonplatten und abgeschlossenen zylindrischen Lufträumen. 
Das streng angesetzte Problem kann mit Hilfe der Ritzschen 
Methode beliebig genau gelöst werden. Weiterhin wird die 
annähernde Berechnung auf einem anderen Wege durchgeführt, 
der sich an die Schwingungen von Massenpunkten anschließt. 
Die Berechnung ergibt die Lage der Eigentöne für beliebige 
Höhen der zylindrischen Lufträume. 


Die vorliegende Arbeit wurde im Physikalischen Institut 
der Universität Breslau auf Anregung und unter Leitung von 
Hrn. Prof. Dr. Waetzmann ausgeführt. Für seine zahlreichen, 
wertvollen Ratschläge und das fördernde Interesse, das er der 
Arbeit entgegenbrachte, bin ich ihm zu größtem Danke ver- 
pflichtet. 


Breslau, Universität, September 1927. 


(Eingegangen 14. September 1927) 
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